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Abstract 



ff^ i The renormalization group (RG) method as a powerful tool for reduction of evo- 

I lution equations is formulated in terms of the notion of invariant manifolds. We 

start with derivation of an exact RG equation which is analogous to the Wilsonian 
, RG equations in statistical physics and quantum field theory. It is clarified that 

I the perturbative RG method constructs invariant manifolds successively as the ini- 

' tial value of evolution equations, thereby the meaning to set to = ^ is naturally 

0^ . understood where is the arbitrary initial time. We show that the integral con- 

stants in the unperturbative solution constitutes natural coordinates of the invariant 
manifold when the linear operator A in the evolution equation is semi-simple, i.e., 
diagonalizable; when A is not semi-simple and has a Jordan cell, a slight modifi- 
r~| ! cation is necessary because the dimension of the invariant manifold is increased by 

the perturbation. The RG equation determines the slow motion of the would-be 
integral constants in the unperturbative solution on the invariant manifold. We 
present the mechanical procedure to construct the perturbative solutions hence the 
initial values with which the RG equation gives meaningful results. The underlying 
structure of the reduction by the RG method as formulated in the present work 
turns out to completely fit to the universal one elucidated by Kuramoto some years 
ago. We indicate that the reduction procedure of evolution equations has a good 
correspondence with the renormalization procedure in quantum field theory; the 
counter part of the universal structure of reduction elucidated by Kuramoto may 
be the Polchinski's theorem for renormalizable field theories. We apply the method 
to interface dynamics such as kink-anti-kink and soliton-soliton interactions in the 
latter of which a linear operator having a Jordan-cell structure appears. 
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1 Introduction 



There is an ever growing interest in the renormahzation groups (RG) fl], in various 
fields of science and mathematical physics since the work of Wilson ^, p. The essence 
of the RG in quantum field theory (QFT) and statistical physics may be stated as follows: 
Let r(0, gr(A),A) be the effective action (or thermodynamical potential) obtained by in- 
tegration of the field variable with the energy scale down to A from infinity or a very large 
cutoff Aq. Here Sf(A) is a collection of the coupling constants including the wave-function 
renormahzation constant defined at the energy scale at A. Then the RG equation may 
be expressed as a simple fact that the effective action as a functional of the field variable 
(j) should be the same, irrespective to how much the integration of the field variable is 
achieved, i.e., 

r(0,g(A),A) = r(0,g(A'),A'). (1.1) 

If we take the limit A' — > A, we have 

c?r(0,^(A),A) 



dA 



0, (1.2) 



which is the Wilson RG equation^, or the fiow equation in the Wegner's terminology 
notice that Eq.(|L2|) is rewritten as 

dg^_dr 

dg dA dA' ^ ' ' 

If the number of the coupling constants is finite, the theory is called renormalizable. In 
this case, the functional space of the theory does not change in the fiow given by the 
variation of A; one may say that the fiow has an invariant manifold. 

A notable aspect of the RG is that the RG equation gives a systematic tool for ob- 
taining the infrared effective theories with fewer degrees of freedom than in the original 
Lagrangian relevant in the high-energy region. This is a kind of reduction of the dy- 
namics. Finding effective degrees of freedom and extract the reduced dynamics of the 
effective variables in fact have constituted and still constitute the core of various fields of 
theoretical physics. In QCD written in terms of the fields of quarks and gluons, the low 
energy effective theories in which gluons are integrated out may be Nambu-Jona-Lasinio 
type lagrangians 0]; see for example If the quarks are further integrated out the ef- 
fective theories are sigma models written solely with meson fieldsQ: The fiow equations, 
variants of the RG equation, may give a foundation on such sigma models as effective 
theories of QCD at low energies p^, [Tl|, |l^; for a review, see |jl3|: In terms of notions in 



the theory of dynamical systems, the functional space represented by a sigma model may 
be an attractive submanifold of QCD. 

Statistical physics may be said to be a collection of theories on how to reduce the dy- 
namics of many-body systems to one with fewer variables, since the work of Boltzmann[Q. 



3 



Bogoliubov showed that BBGKY( Bogohubov-Born-Green-Kirkwood-Yvon) hierarchy can 
be reduced to Bohzmann equation with a single-particle distribution function for dilute 
gas systems. ||15|| As indicated by Kuramoto||16||, Bogoliubov seems to have claimed that 



the dilute-gas dynamics has an attractive manifold spanned by one-particle distribution 
function. Boltzmann equation in turn can be further reduced to the hydrodynamic equa- 
tion (Navier-Stokes equation) by a perturbation theory like Chapman-Enskog method [p!7[]. 
Recent development of the theories of pattern formation with dissipative structures gives 
a good example how to reduce complicated ordinary and partial differential equations to 
simple equations with slow variables, such as Landau-Stuart equation, the time-dependent 
Ginzburg-Landau equation and so on. 



Some years ago, it was shown by an Illinois group |21|, ^ and Bricmont and 



Kupiainen ||24|| that the RG equations can be used for a global and asymptotic analysis of 
ordinary and partial differential equations, hence giving a reduction of evolution equations 
of some types. A unique feature of the Illinois group's method is to start with the naive 
perturbative expansion and allow secular terms to appear; the secular terms correspond 
to the logarithmically divergent terms in QFT. Then introducing an intermediate time 
T and rewriting perturbative solutions in terms of renormalization constants reminiscent 
of those appearing in the perturbative renormalization theory in QFT, Gell-Mann-Low 
type RG equation is applied to obtain the evolution equations for the renormalization 
constants which are functions of r. Finally, equating r with the time t appearing in 
the original perturbative solution, they obtained global solutions of differential equations. 
Bricmont and Kupiainen p4[ applied a scaling transformation (block transformation) to 



obtain asymptotic behavior of nonlinear diffusion equations in a rigorous manner. 

Subsequently, one of the present authors (T.K.) formulated the Illinois group's method 
in terms of the classical theory of envelopes |2^, |26[: He indicated that the RG equation 
can be interpreted as the basic equation for constructing envelopes of a family of curves (or 
surfaces for partial differential equations). He also developed a short-cut prescription for 
the renormalization procedure without introducing an intermediate time r but utilizing 
an arbitrary initial time to. In latest papers |2^, it was shown that the RG method can 



be well formulated as the method dealing with the initial values at arbitrary t = Iq] the 
initial values at t = to are determined so that the unperturbative solutions which are 
valid only locally around t = to are continued smoothly; this procedure is nothing but to 
construct the envelope of the perturbative solutions. He emphasized that the RG method 
is a powerful tool for reduction of evolution equations and demonstrated it by applying 
the method to obtain so-called the amplitude equations for systems of equations. He 
also suggested that if the phase equations describing slow motions in the system where a 
continuous symmetry is broken, which are classical counter part of the Nambu-Goldstone 
modes in QFT, the RG method should be able to derive the phase equations. 

The RG method developed by the Illinois group has been applied by many authors 



to quite a wide class of problems successfully. To mention some of them; Graham p8 
derived a rotationally invariant amplitude equation appearing in the problem of pattern 
formation. Some kinds of phase equations were also derived; Oono p9[ discussed the 
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interface dynamics relevant in the spinodal decomposition, Sasa derived a diffusion type 
phase equation, and Maruo et al ||31| derived Kuramoto-Sivashinsky equation[0. The 
method was apphed to analyze asymptotic behavior of the non-linear equations appear- 
ing in cosmology p3|, EM. Boyanovsky and de Vega also used the method to derive an 



anomalous transport coefficient relevant in the non-equilibrium states in the early 
universe and QGP (quark-gluon plasma). Before that, the Boltzmann equation had been 
derived as an RG equation [^. The RG method was also shown to be a powerful tool to 
resum divergent perturbation series appearing in problems of quantum mechanics P7|, 
Tzenov applied the method to obtain global solutions appearing accelerator physics ||38[ . 
Possible relation between renormalizability and integrability of Hamilton systems was 
discussed by Yamaguchi and Nambu||39[|. The method was proved to be applicable to 
discrete systems too. 



Although such extensive applications have been made, only few works are known 
which attempt to reveal the underlying reasons why the RG method works to some kinds 
of equations but not to others[^]. Nevertheless, it was indicated in ^] that the 
RG method by the Illinois group works when the unperturbed solutions are neutrally 
stable solutions that are stationary (constant in time or stationary oscillation) hence do 
not decay nor blow up with time. A decade ago, Kuramoto revealed in an excellent 
paper [0 the universal underlying structure of all the existing perturbative methods for 
reduction of evolution equations; he noticed that when a reduction of evolution equation 
is possible, the unperturbed equation admits neutrally stable solutions, and succeeded 
in describing the reduction of dynamics in geometrical terms, i.e., attractive manifolds 
or invariant manifolds]^. Although his actual presentation of the theory was based on 
the reductive perturbation theory and involves some ansatz on forms of the solutions, 
he emphasized that the universal structure which he revealed should not be dependent 
on the perturbation methods one employs. In the present paper, focusing on the aspect 
of the RG method as a powerful tool for reducing evolution equations, we shall present 
a comprehensive formulation of the perturbative RG method in terms of the notion of 
invariant manifolds, guided by this Kuramoto's work: One will see that his ansatz are 
derived naturally in the RG method. It may mean that his formulation of reduction 
of evolution equations, which is actually a natural extension of the asymptotic method 
by Krylov, Bogoliubov and Mitropolski for non-linear oscillators HTI 



is an RG theory 



although the term RG is not usedp6| . 



We start with adapting the exact RG equations of Wilson type (flow equations) in 
quantum field theory || to differential equations (evolution equations); one may rec- 
ognize that the RG method described in terms of envelopes is best formulated in the 
framework of Wilson RG. Then confining ourselves to cases where a perturbative treat- 
ment is possible, we shall show that an invariant manifold exists when the RG method 
works, and that the RG method is a method to construct the invariant manifold and the 
reduced dynamics on it in a mechanical way. In this method, the initial values of the so- 
lution at an arbitrary time t = in the successive perturbation orders are determined so 
that the initial value make an invariant manifold, thereby the condition to set Iq = t will 
be naturally emerged. We emphasize that the present formulation give a nice foundation 
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for the prescriptions adopted in pSI, EBI, P? . 



The following should be mentioned here: (1) The relevance of Wilson type RG equation 
to their method was noted by Chen et al[^ and Pashkov and Oono p6| , although their 
formulation was totally based on the Gell-Mann-Low type perturbative RG method. (2) 
Shirkov ^ clarified that the RG equation concerns with the initial values and emphasized 
the Lie group structure of the RG method ||4^; he extracted the notion of functional self- 



similarity (FSS) as the essence of the exact RG. He claims that the Wilson RG is an 
approximation to the Bogoliubov RG which is exact 0. 

Once the underlying structure of the reduction of dynamics given by the RG method 
has become clear, one will recognize that the method could be applied to problems which 
have not been treated in the RG method so far although other methods are applied 
to them. One will also see that the RG method is simpler to apply to them than the 
previous methods. As such a problem, we take the problem to extract the interface 
interactions]^, ^ |45| , which are typical examples to be treated by the method of 



phase equations. 



In §2, we formulate the RG method as a method of reduction of dynamics, starting 
from the non-perturbative flow equation (RG equation). Then on the basis of the per- 
turbation theory, we show, guided by the presentation given in |]16|, the way how the 



invariant manifold which is supposed to exist to the evolution equation under considera- 
tion can be constructed in our method. We shall remark that the existence of an invariant 
manifold corresponds to the notion of the renormalizability in quantum field theory. In 
§3, some simple but typical examples including the Takens equation (091 are worked out 



to demonstrate how the RG method construct invariant manifolds and give the reduced 
dynamics on them. In §4, we show how to deal with generic systems which involve a linear 
operator A having zero-eigenvalue where A may or may not be diagonalizable; when A 
is not diagonalizable, the eigenvalues are degenerate and A is equivalent with a matrix 

1 



having a Jordan cell 







As examples with a Jordan cell, we shall show that the 



RG equation gives the normal forms of the two-dimensional equations including the 
Takens and the Bogdanov equation and deal with an extended version of Takens 
equation with three-degrees of freedom. In §5, we apply the method to some problems 
such as the unstable motion in the Lotka-Volterra system and the Hopf bifurcation in the 
Although the examples treated in 



Brusselator 52 



^2-5 are simple ordinary differential 
equations, we believe that these examples should be instructive also for experts of the 
RG's or flow equations in quantum field theory and/or statistical physics. In §6, we also 
apply the method to extract the interface dynamics of a kink-anti-kink interaction in the 
TDGL equation in one-dimension and the soliton-soliton interaction in the KdV equation. 
The final section is devoted to a brief summary and concluding remarks. In Appendix A, 
we summarize rules in a scheme of an operator method for constructing special solutions 
suitable for the RG method. In Appendix B, we present an elementary method different 
from the RG method to derive the approximate solution for the double-well potential dis- 
cussed in subsection 3.3. In Appendix C, we show that the period of the Lotka-Volterra 
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system around the non-trivial fixed point obtained previously[^ in the RG method 
coincides with that extracted by Frame in a quite different way. 
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2 Reduction of evolution equations with the RG method 



In this section, we shall formulate the RG method in terms of the notion of invariant 
manifold. 



2.1 The RG equation as the flow equation for initial values 



2.1.1 Non-perturbative RG equation 



Let us take the following n-dimensional dynamical system 

dX 



dt 



F(X,t), (2.1) 



where n may be infinity. We remark that F may depend on t explicitly. We suppose that 
the equation is solved up to an arbitrary time t = to from an initial time, say, at t = 0, to 
give X(t) = W(t), and then we are trying to solve the equation with the initial condition 

at t = Vto, 

X{t = to) = W{to), (2.2) 

with Wito) being unspecified yet. In fact, W(t) as a function of t is a solution of ( p.l|) by 
definition. The solution may be written as X(t;tQ, W(to)). We stress that the solution 
needs not be given by a perturbation method; if the solution is given non-perturbatively, 
the resultant equations remain non-perturbative ones. 

Now, making use of X{t]tQ,W{to)) thus obtained, we determine W{to) based on a 
simple fact of differential equations. When the initial point is shifted to tg, the resultant 
solution should be the same, i.e., 

X{t- to, W{to)) = X{t- 1'„ W{t'o)). (2.3) 



Taking the limit tg to, we havep7 



dX dX dX dW , 

= I = f2 4) 

dto dto dW dto ' 

This equation gives the evolution equation or the flow equation of the initial value W(to). 
This has the same form as that of the renormalization group (RG) equation in quantum 



field theory, hence the name of the RG method. We emphasize that the equation ( |2.4| ) 
is exact; we have not used any argument based on perturbation theories. This equation 
corresponds to the non-perturbative RG equations (fiow equations) by Wilson|Q, Wegner- 
Houghton^] and so on in quantum field theory and statistical physics|]l3|: The reader 



should have recognized that to corresponds to the logarithm of the energy scale to which 
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the integration is performed in the quantum field theory. The quantities corresponding 
to the couphng constants will be found to be the integration constants. One will also 
recognize that the existence of an invariant manifold of a dynamical system corresponds 
to the renormalizability. We also notice that one needs not to equate to with t at this 
stage. 



2.1.2 Perturbative RG equation 



So far, only the perturbative expansion method is available to make X{t]to,W{tQ)). In 
this case, X(t;to,Wito)) and X{t]t'Q,W(tQ)) may be vahd only for t ~ to and t ~ tg. 
This condition is naturally satisfied if we restrict that < t < t'^ (or I'q < t < to) 
because the limit tg — * to is to be taken. Thus when a perturbative expansion is used for 
constructing X{t; to, W{to)), we demand a more restrictive RG equation as given by||27 



dX 



dto 



dX 



to=t 



dtn 



dX dW 



to=t 



dW dto 



0. 



to=t 



(2.5) 



Notice that the demand to set to = t has naturally emerged. 



We remark that a geometrical interpretation of this equation has been given on the 
basis of the classical theory of envelopes |2^, When to is varied, X(t;to, W(to)) 
gives a family of curves with to being a parameter characterizing curves. Then Eq. (p.5|) is 



the condition to construct the envelope of the family of curves which are valid only locally 
around t ~ tg. The envelope is given by X(t;to = t) = W{t), i.e, the initial value. It is 
noteworthy that W{t) satisfies the original equation (|2.1| ) in a global domain up to the 
order with which X{t; to) satisfies around t ~ to [p5| , p6| . In fact, one can see that 



dW 

IT 



dX{t-to) 



dt 

dX{t;to) 



dt 



dX{t-to) 



to=t 



to=t 



dtn 



to=t 



(2.6) 



on account of (0| 



2.2 Invariant manifolds and renormalizability 



In this section, we follow ||TB| for notations. If the theory has an invariant manifold M 
with the dimension less than n, we may have supposed that the initial point is on the 
manifold. Let the invariant manifold M is represented by the coordinate s. The reduced 
dynamics of Eq. ( |2.1| ) on M may be given in terms of a vector field G by 

I^GW. (2.T) 
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and the manifold M is represented by 

X = R{s). (2.8) 

Our task is to obtain the vector field G and the representation of the manifold R in a 
perturbation method. We consider a situation where the vector field F is composed of 
an unperturbed part Fq and the perturbative one P, i.e., 

F = Fo{X) + e-P{X,t). (2.9) 

Here notice that Fq has no explicit t-dependence, while P{X,t) does. We assume that 
the unperturbed problem is solved and an attractive invariant manifold Mq is easily found. 



Now we try to solve Eqs.(2J.) and ( p.9|) by a perturbation theory with the initial 
condition 

X{to) = W{to), (2.10) 

at t = to- The decisive point of our method is to assume that 

W{to) = R{s{to)), (2.11) 

that is, the initial point is supposed to be on the invariant manifold M to be determined. 
Now we apply the perturbation theory, expanding 

X{t; to, W{to)) = Xo + Yl ^0, W{to)). (2.12) 

Here we have made it explicit that X is dependent on the initial condition. We should 
also expand the initial value, 

W{to) = Woito) + pito), (2.13) 

with 

oo 

P(^o) = Ee"W^-(^o). (2.14) 

n=l 



Now the unperturbed equation reads 

dXo 



FoiXo). (2.15) 



dt 

As promised, we suppose that an attractive manifold is found for this equation as 

Xoit) = R{s{t;C)), (2.16) 
where C is the integral constant vector with dimC = m and may depend on to. 
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Here comes an important point of our method; we identify that 



s{to) = C{to), (2.17) 

which gives a natural parameterization of the manifold Mq. This is a simple but a sig- 
nificant observation; notice that we need not to give any ansatz to the representation of 
the manifold because we only have to solve the unperturbed equation and the integral 
constants are trivially obtained. 

The deformation of the manifold p will be determined perturbatively on the two 



principles, i.e. 16 



1. the function p should be independent of Wo, and 

2. the resultant dynamics should be as simple as possible because we are interested to 
reduce the dynamics to a simpler one. 



The choice of the p is intimately related to that of the forms of the perturbative 
solutions. For example, the first order equation reads 

^ = n(^o)Xi + P(Xo). (2.18) 

The solution to this inhomogeneous equation is composed of a sum of the general solution 
of the homogeneous equation and the special solution of the inhomogeneous equation. 
If the unperturbed solution Xo(t) has a part of neutrally stable solution, there appear 
secular terms in the special solution as well as genuinely independent functions. It turns 
out that the secular terms can be utilized to renormalize out the homogeneous solutions 



at t = to[il; |26[; this is a kind of renormalization conditions. Then the shift of the initial 



value is now determined as 



W,{to) = X^{t = to). (2.19) 

Notice that the initial value is determined after solving the equation, hence the functional 
form of it as a function of to is explicitly given. Of course, one can make Wi(to) = 
by further adding unperturbed solutions; this prescription, however, will generally give 
a more complicated dynamics. It means that the choice of p has ambiguities, but the 
demand to obtain the simplest dynamics give it uniquely. In §3 and 4, we shall give a 
mathematical and mechanical procedure to select the initial values in accordance with 
the above rules using some examples. 



We can repeat the procedure to any order of the perturbation. We remark that by 
this procedure the modification of the initial value p{to) and hence the total initial value 
W(to) are given solely in terms of C(to); 

W{to) = Wo[C]+p[C]. (2.20) 
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Now the dynamics of C(t) is determined by the RG equation Eg. ; 



dX 

dtn 



to=t 



dX 



to=t 



dX dC 
"dC ' dt^ 



0, 



to=t 



(2.21) 



which is an evolution equation of C{t). Then the manifold M (more precisely, the trajec- 
tory on it) is represented as 



X = Wit) = Wo[C{t)] + p[C]. 



(2.22) 



Eq.'s ( [^.21| ) and ( |2.22| ) are essentially the basic postulates in Kuramoto's theory of reduc- 
tion of evolution equations |jl6|. Thus we have derived the Kuramoto's basic equations in 
the RG method; in other words, Kuramoto's theory is actually the RG theory for reduc- 
tion of evolution equations. Eq.( |2.22| ) means that the dynamics is renormalized to that 
for C. One can now see a correspondence between the renormalizability of the theory in 
quantum field theory|p and the existence of a finite dimensional invariant manifold in the 
theory of dynamical systems: C = (Ci, C2, ■ ■ ■ , Cm) correspond to the collection of renor- 
malizable coupling constants and p to unrenormalizable operators; one may say that the 
fact that p can be represented solely with C is analogue of the Polchinski theorem]^, Q 
in quantum field theory. We remark also that Eq.( |2.22| ) justifies the slaving principle by 
Haken||]. 



A comment is in order; When the unperturbed system is given by a linear operator 
with a Jordan cell, there will be a slight modification of the above scenario due to a 
technical complexity; see an example in the next section and the general argument given 
in 54. 



3 Simple examples 



In this section, we consider four simple equations to show our formulation of the RG 
method at work as a tool for reduction of evolution equations. We shall show how the 
initial values at t = Vto in higher orders are determined by the two principles that terms 
proportional to the unperturbed solution are suppressed and that possible fast motions 
disappear. The resultant forms of the unperturbed solutions composed of the secular 
terms which vanish at t = and the functions independent of the unperturbed solution. 
The final forms turn out to be the same as those given in the scheme adopted in ||25|, ^ p^ , 
where unstable manifolds and cases with a Jordan cell were not dealt with though. The 
resultant initial value represent the invariant manifold, and the RG equation gives the 
slow dynamics on the manifold. 
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3.1 A simple model with a reduction 



We first consider the following simplest equation with a reduction||16||; 



dx/dt = ef{x,y), dy/dt = —y + g{x). (3.1) 

Writing u{x, y) = *(a;, y) and we expanding u as u = uq + eui + ■ ■ ■ , with Un{x, y) = 
^{xn,yn) {n = 0, 1, ...). We solve the equation with the initial value W(to) at t = to. We 
suppose that W{to) is on an attractive manifold M. 

The unperturbed equation reads 

Xo = 0, i/o = -yo + 9{xo), (3.2) 

the solution to which is readily obtained as xo(t) = const. = Co, yo(t) = g{Co) + Cie~*, 
with Ci {i = 0, 1) being the integral constants; we make it explicit that Ci may depend 
on the initial time to as Ci = Cj(to)- Since we are interested in the asymptotic behavior 
as t — >• cxo, we take the stationary solution putting Ci = 0; 

xo{t) = const. = Co, yo{t) = g{Co), (3.3) 

accordingly, 

Woito)= *(Co,^?(Co)), (3.4) 
which suggests that the unperturbed invariant manifold Mq is given by 

y = g{x), (3.5) 
although the time dependence x{t) is not yet known. 

The first order equation reads 

(3,'A)H = {f^^''^«A (3,6) 



with 

A-( 

We notice that 



AU, = 0, AU2 = (-1)1/2; C7i = , = (^^ • (3.8) 
Then the special solution with the initial value W^i(to) is obtained as follows: 



h] = e(*-*«)-^i¥i(to)+ r 

\yi J J to 



dse^'-'^^l 



(f{xo,yo)\ 

e(*-*°)-^[iyi(to) + fg'U^] + {{t - to)/C/i - fg'U2}. (3.9) 
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Here we have used the relation that 0) = Ui — g'U2- For the solution to describe a 
slow motion, the first term should vanish; thus we are naturally led to the choice of the 
initial value as 

W^{to) = -fg'U2, (3.10) 

to kill the fast motion. Thus we have also 

xi{t- to) = /(Co, giCoW - to), yiit- to) = /(Co, g{Co))g\Co){{t - to) - 1}. (3.11) 

We notice that the solution and the initial value satisfy the rules given in the preceding 
subsection; Wiito) — p(to) is independent of W^o(^o)- 

Up to this order, 

„(f.,^_{<t-^to)\ _ ( Co{to) + ef{Co,g{Co)){t-to) \ . . 

^ ' [ yit; to) ) [ g{Co) - eg'{Co)f{Co, g{Co)){l + to - t) J ' ^ ^ ^ 



with 



W{to) = f^^JJ^I) ^ Wo{to) + p{to) = f _ , _ , ,) . (3.13) 

\yv-G)) \9{Co) - eg'{Co)f{Co, g[Co))) 



Now the RG equation Eq.(|2.5|) applied to Eq.( |3.12|) gives the evolution equation for 
Co(t); 



dCo 
~dt 



tf{Co,g{Co)). 



(3.14) 



Since x{t) = Co(t), Eq. (|3.14|) gives the reduced dynamics, and the slow manifold is given 
by u{t) = W{t), or in terms of the components. 



x{t) = Co(t), y{t) = g{Co) - eg'{Co)f{Co,g{Co))). 
One sees that the attractive manifold is given by 

y{x) = g{x) - eg'{x)f{x,g{x)), 



(3.15) 



(3.16) 



and the original two-dimensional evolution equation has been reduced to the one-dimensional 
equation Eq.( |3.14|) . 



In short, to obtain the invariant manifold and the slow dynamics on it, we have started 
with a neutrally stable solution with an integral constant, which constitutes a natural 
representation of the invariant manifold. The constant moves slowly by the perturbation 
and the whole dynamics is given through this slow variable. The evolution equation of the 
slow dynamics is given by the RG equation and the trajectory on the invariant manifold 
is given as the initial value in the RG method uniquely. The initial value in the higher 
order has been determined so that the fast motion disappears. 

Generic systems which have a linear matrix having zero eigenvalue will be extensively 
analyzed in §4. 
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3.2 When the unperturbed solution is oscillatory 



In the present subsection, we shall examine a case where the unperturbed solution has 
another type of neutral stability, i.e., is oscillatory . We shall treat a simplest equation 
of a damped oscillator with the geometrical terms putting an emphasis on the fact that 



the RG method concerns with the initial value. In a recent monograph, Nishiurap? 
used this example to indicate that a careful identification of the initial values in the 
higher order terms is needed to construct the proper perturbative solutions for obtaining 
a meaningful result in the RG method, though he failed to give general principles for "a 
careful identification". One will see that the general principles given in the foreword in 
this section determine the initial values uniquely, hence Nishiura's concern is resolved. 
One will also see that the initial values thus obtained are of the same forms as obtained 
in 
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The equation we deal with is 

x + ex + x = 0, (3.17) 

where < e < 1. This system does not exhibit a decrease of the degrees of freedom, 
nevertheless the dynamics is reduced to a set of simpler equations for the amplitude and 
the phase, separately by the RG equation. 



With the definition u = ^{x,y), y = x, Eq.( p.l7| )is converted to 



with 



(dt-A)u = -e{^], (3.18) 



^ = ( ° 1 ) • ^3-^9) 



We expand the dependent variable and the initial value in Taylor series as u{t) = 

uq + eui + e^U2 H and Wito) = Wq + eWi + e^W2 H , with Ui = \xi, yi). The 

lowest order solution reads with the (complex) integral constant C(to); 

Uo{t-M = Ce^'U^ + C%-^'U. ^ (3.20) 

where U± = *(l,±i); AU± = ±iU±. This is a neutrally stable solution. The initial 
value in this order reads, 

Woito) = z{to)U+ + c.c., (3.21) 

with 

zito) = C(to)e^*°. (3.22) 
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Here c.c. denotes the complex conjugate. 
The first order equation reads 



1 



{dt-A)u, = -yo{t)-{U+-UJ, 



(3.23) 



with the initial condition ui{to, to) = Wi{to) which is not yet known but to be determined. 
The equation is readily solved as follows; 

u,{t;to) = e(*"*°)^Wi(to)-T^ rrfse(*-^)^o(3)(C/+-C7_), 

2i J to 



WAto) {C*e~**«C/+ - Ce**«C/_} 

1 



{{t - to)Ce''U+ - —Ce^'U^} + c.c. 



(3.24) 



which leads to the natural choice of the initial value 

Wiito) = ^^{z*{to)U+-z{to)U^} = p,[z,z*], (3.25) 

because otherwise the first term of ( p.24|) gives rise to terms which could be renormal- 
ized away into the unperturbed solution with a redefinition of C. We emphasize that a 
renormalization procedure enters here. Hence 

'xi(t;to)\ 



ui{t; to) = -liit - to)C{to)e''U+ - lc(to)e'*C/_} + c.c. = . , , 
/ zt \yi[t,to) 



(3.26) 



Similarly, the second order solution is given by 

1 /•* 



li ho 



^{t-to)A 



iC 



Vr2(to)-{-^e**°C/- + c.c.} 



+— [{z(t - to) + (t - to)}C/+ + {t - to + V2}C/-]e' + c.c. (3.27) 



Thus we are led to the choice of the initial value 

^2(^0) = -^z{to)U- + c.c. = p^\z, z* 
lb 

because of the same reason as in the first-order case. Hence 



tt2(t; to 



iC 



[{i{t - to f + (t - to)}t/+ + {t - to + z/2}t/_]e^* + c.c. 



(3.28) 



(3.29) 



Collecting all the terms up to the second order, we have 
w(t;to) = Ce''U^ 



^{(t-to)C/+ + ^C/_}e^ 



+ ^[{z(t - to)' + (t - to)}C/+ + {t - to + t/2}U.]e'' + c.c, (3.30) 
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with the initial value 



W{to) = z{to){U+ - t-U. - -t/_} + c.c. (3.31) 

4 lb 



The RG equation Eq.( |2.2lD gives 

(7 + (e/2 + zeV8)C = 0. (3.32) 

Parameterizing C as C = {A/2) ■ exp{i6), we have the equations governing the amplitude 
and the phase, respectively, 

A + e/2-A = 0, 9 = -e'^/8, (3.33) 

which yields A{t) = Aqc"''/^'* and 6{t) = — e^t/S + 6*0 with Aq and Oq being constants. 
Since n(t) = W{t), we have the final solution to the damped oscillator as 

x{t) = Aoe-'/^-'{{l - — ) cos{ut + 9o) + - sin(u;t + ^o)}, (3.34) 

lo 4 

with u = 1 — being the angular velocity. The above expression is slightly different 
from that given in where the equation is treated as a scalar equation; a redefinition 



of the constants Aq and 9q transforms the solutions to each other in this order. We see 
that although the number of the dimension of the equation is not changed, the dynamics 
is reduced to simpler equations ( p.33| ) for the amplitude and the phase. The initial values 



in the higher order equations have been determined so that terms proportional to the 
unperturbed solution do not appear; such higher order terms have been "renormalized 
away" by a redefinition of the integral constant in the unperturbed solution. 

A comment is in order: The above solution could be more efficiently obtained by using 
the operator method (implicitly) adopted in |2^, ^ and fully accounted in Appendix A 



of the present article. For instance, the first order solution in the operator method reads 
Hi(t; to) = -^idt - A)-\Ce'' - C*e-**)(L/+ - t/_), (3.35) 

which is readily evaluated with the use of (A. 9) and (A. 10). The second order solution 
is also readily obtained with the use of the formulae given in Appendix A. 



3.3 Unstable motion in the double-well potential 

Next, we show that unstable manifolds are also constructed by the present method, using 
the double-well potential in mechanics; 

X = X — ex^. (3.36) 



17 



Let us obtain the unstable motion around the origin by applying the RG method. This 
example is a peculiar case where the unperturbed solution is not neutrally stable but 
composed of a blowing and decaying function with the same exponent. We shall treat 
another example of this type in §5.2. 



Putting X = y and defining u = ^{x,y), we have a system of equation 



-A]u = e[ A. (3.37) 



We solve the equation around an arbitrary t = with the initial value VF(to)- The 
solution is written as u = u{t; to, W{to)). We expand u = Uq + eui + ■ ■ ■ . Accordingly 
the initial value Wito) which is to be determined self-consistently is also expanded as 
W ^Wq + eWi + ■■■. 

The lowest order solution reads 

Uo{t; to) = C+{to)e'U+ + C_(io)e-*C/_, (3.38) 

where we have suppressed the W-dependence of Uq and U± — *(1,±1) are the eigen- 
vectors of A belonging to the eigenvalues ±1, respectively. The initial value Wo(io) 
accordingly reads 

Woito) = Wo(*o; *o) = C+(to)e*°C/+ + C_(io)e-*°C/-, (3.39) 
which will imply that the trajectory is in a hyperbolic curve 

Mo = {u = {x, y)\{x + y) {x - y) ^ const} (3.40) 
with the asymptotic lines y — ±x. 

The first order equation reads 

(I - ^) Hi = -\{C+{toy + C^{h)e-'f{U^ - C/_), (3.41) 

which is solved with the initial value Wiito) formally as 
u^{t-to) = e(*-*°)^Wi(to) 

_1 f dse^'-^^\C+{t,)e^ + C_{t,)e'^ f{U+ - L/_), 

= e(*-*°)^[Wi(io) + to)U+ + f4to; to)U_)] 

-^{U{t;to)U+ + f4t;to)U_), (3.42) 



with 



/+(t;to) = 3ClC4t-to)e'-^C^Cle-' + ^Cy-^Cie-'\ (3.43) 
f4t;to) = -3C+C'_{t-to)e-'-^ClC^e' + ^Cie-''-^Cle''. (3.44) 
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Thus the first-order initial value Wi(to) is chosen to be 

1 



;(/+(to;to)C/+ + /-(to;to)C/-), 



(3.45) 



because otherwise the first term of ( p.42|) gives rise to terms proportional to the unper- 
turbed solution which should be "renormalized away" with the redefinition of C±. We 
notice that Wiito) is a function of A±{to) = C±(to)e^*°; 



W,{to) = pAA+,A_]. 



(3.46) 



Applying the RG equation to u(t; to) = Uo{t; to) + eui(t; to), one obtains 



du 

dto 



= C+e'U+ + C.e^'U^ + 3-{C2 C_e*C/+ - C+C!e-*C/+} = 0, (3.47) 
to=t 2 



which leads to 



3e 
2 



3e 
~2 



(3.48) 



Noting that C+(t)C_(t) = const = c+c_, one obtains the solution to the RG equation as 
follows; 



C+{t) = c+e~^''+'-'/\ C_{t) = c_e=^^^+^-*/2_ 
Thus one finds that the global solution is given by 

uit) = W{t) = Wo[A+,A_]+ep,[A+,A_], 
= A+{t)U+ + A^{t)U- 

-l[{2Al{t) - QA4t)Al{t) - Al{t)}U^ 
+{-Alit) - 6Al{t)A4t) + 2Al{t)}U^ 



(3.49) 



(3.50) 



where A±(t) = c±exp{±at} with a = 1 — 3ec+c_/2. In this case, the unstable manifold 
is represented solely with A± and the dynamics on the manifold is given through the 
evolution of these variables. 

To compare the result with the exact solution given in terms of an elliptic function, 
we first write down the first component x{t) of u{t) given in ( |3.50| ); 



x{t) = (1 + -e/3)(c+e"* + c_e-"*) - -{^ e^"* + c^e-^"*}, 
4 8 

with P = c^c_. Writing (1 -|- je[3)c± as c±, we have 



(3.51) 



x{t) = c+e"* + c_e 



-at 



(3.52) 
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still with a = 1 — 3ec+c_/2 up to O(e^). In fact, if one solved the equation in the scalar 
form without converting the equation to the vector one and applied the operator method 



adopted in ||2^, ^ and explained in Appendix A of the present paper, one would have 
directly reached the form given in ( p.52| ). We stress that the method adopted in [^S], ^ 
is more convenient in practice than that presented here. 

To make the comparison with the exact solution easier, it is found convenient to 
consider the case with the initial condition x(0) = 0. This implies that c+ = — c_. 
Further putting c+(l + 3e/3/8)/2 = C, one has 

x{t) = C sinh at - -C^ sinh^ at, (3.53) 
8 

with a = l + feC^ q(^^2-^_ 

Now the exact solution with the initial condition x{0) = reads 

x{t) = h cn{a{e)t + K{k),k), 

= -hk'F{a{e)t,k), (3.54) 



where F{a{e)t,k) = sn{a{e)t,k)/ dn{a{e)t,k) with a{e) = yje/2k/h, cn(t, k), sn(t, k) and 
dn(t, k) are Jacobi's elliptic functions with modulus fc, K{k) the complete elliptic integral 
of the first kind and k' = yl — W the complementary modulus. The constants k and h 
are functions of e and the energy E of the system, which is assumed to be positive here; 



and a{e) ^ 1 + eE. Notice that k approaches 1 as e goes to 0. If one expands F(Q;(e)t, k) 
w.r.t. k around k = 1 with a{e) fixed, one has 

F{a{e)t, A;) ~ (1 + ^E) sinh(l + -eE)t - ^E sinh^(l + -eE)t, (3.56) 

up to O(e^). Notice that the expansion of the elliptic functions at /c = 1 is subtle; we 
have made the manipulation as follows; 

^nhu-'-Eucoshu = sinhu-'-Eu^^, 
4 A du ' 

sinh(l - ^E)u. (3.57) 



Then identifying C = —hk'{l + ^E) = —V2E{1 + |e-E), one reproduces the above result 
( |3.53|) . In Appendix B, an elementary method is presented to derive the approximate 
formula obtained above. 
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3.4 An example with a Jordan cell; Takens equation 



As examples which involve a linear operator with a Jordan cell, we take the Takens 
equation |^ in the present subsections. 



The Takens equation is given by 

X = y + ax^ , y = bx^. (3.58) 

Since we are interested in a slow motion in the vicinity of the origin, we make a scale 
transformation; 



X = e"X, y = e''Y, (3.59) 

where e is a small parameter. To make the equation to be balanced, one finds that a = (3\ 
we choose a = (3 = 1 for simplicity. Then we end up with 

X = Y + eaX^ Y = ebX\ (3.60) 



Expanding as X = Xq + eXi + ■ ■ ■ , and F = Yq + e^i + ■ ■ ■, we first solve the equation 

around t ~ to with the initial value W{to) = \X{to), Y{to)) = Wo{to) + eWi{to) H , 

where to is arbitrary.i The equations in the first few orders read 

Xo = Fo, Yo = 0, (3.61) 
Xi = Y^ + aXl Yi = bXl (3.62) 

and so on. 

We take the stationary solution as the lowest order one to describe a slow motion on 
an invariant manifold, in accordance with the previous treatment; namely, 

Xo{t- to) = const = Co(to), Yo{t- to) = 0, (3.63) 

Co(to) is an integral constant. Accordingly, Wo(to) = *(Co(to))0), i.e., the unperturbed 
manifold Mq is the X axis. 

Then the first order equation is solved successively from Yi to Xi to yield 
Fi(t;to) = bC^oito){t-to) + Ci{to), 

X,{t-to) = \cl{to){t-tof + {C,{to) + aCl{to)){t-to), (3.64) 

where Ci(to) is another integral constant. Accordingly, Wi{to) = *(0,Ci(to)); namely, 
the modification of the invariant manifold is given in the Y direction; 

Ml = {{X,Y)\iX,Y) = (Co,Ci)}. (3.65) 
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Up to this order, we have 



X{t;to) = Co{to) + e{^C^ito){t-tof + {C,{to) + aC',{to)){t-to)}, (3.66) 

y(t;to) = e{6Co'(to)(t-to) + Ci(to)}, (3.67) 
withW{to)= \Co{to),eC,{to)). 



Now applying the RG equation to X{t;to) and Y{t;to) thus obtained, we have 

Co = e{Ci + aCo{t)), C, = bCl 



(3.68) 



The trajectory is given by 

X{t) = X{t;t) = Co{t), Y{t) = Y{t;t) = eCi(i), 



(3.69) 



which shows that the original Takens equation is reproduced. This means that Takens 
equation is "irreducible" and can not be reduced to a simpler equation. 

A few comments are in order: In this case, the invariant manifold is represented with 
two variables in accordance with the dimension of the Jordan cell although the dimension 
of the unperturbed solution is one; namely, the dimension of the invariant manifold is 
increased from that of the unperturbed manifold Mq. The amplitude of the trajectory in 
the second direction, eCi(t) is small compared with the amplitude in the first direction 
Co, while the time dependence of the second variable is large in comparison with the first 
variable. 



4 Generic systems with the Hnear operator having 
zero eigenvalues 



In this section, we shall examine invariant manifolds and slow motions given by generic 
systems which have a linear operator having zero eigenvalues. We shall show how uniquely 
the initial values are chosen by using a simple formula for the special solutions to differ- 
ential equations as in the previous section; the initial values are determined successively 
so that terms which give fast motions and those proportional to the unperturbed solution 
do not appear. We call these unwanted terms "dangerous ones". We shall also show a 
necessary condition on the type of equations for the RG method to be applicable, which 
condition is relevant when the linear operator has a Jordan cell. 

We treat the following rather generic vector equations in this section: 



where dtU = du/dt, A is a linear operator, F a nonlinear function of u and e is a small 
parameter (|e| < 1). We assume that A has multiply degenerated zero eigenvalues and 
other eigenvalues of A have a negative real part. 




(4.1) 
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We are interested in constructing the attractive manifold M at t — ^ oo and the reduced 
dynamics on it. We try to construct solve the problem in the perturbation theory by 
expanding u as 

u{t; to) = Uo{t; to) + eui{t; to) + e^U2it; to) + ■ ■ ■ , (4.2) 

with the initial value W{to) at an arbitrary time to- We suppose that the equation has 
been solved up to t = to and the solution has the value W{to) at to- Actually, the 
initial value must be determined by the perturbative solution self-consistently; indeed, 
u{t) = W(t) is the solution to (|4.1|) in the global domain. Therefore it should be also 
expanded as follows; 

Wito) = Woito) + eW,ito) + e'W2{to) + ---, 

= Woito) + pito), (4.3) 

where p{to) is supposed to be an independent function of Wo- They are not yet known 
at present but will be determined so that the perturbative expansion becomes valid. One 
of the main purposes in this section is how sensibly the initial values can be determined 
order by order. 

The equations in the first few orders read 



where 

n 

{F'{uo)u,), = {d{F'iuo))^/diuo)j} (mi)„ (4.7) 



if u is an n-dimensional vector. 



{dt~A)uo = 0, 
idt-A)ui = F(no), 
idt-A)u2 = F'{uo)u,, 



(4.4) 
(4.5) 
(4.6) 



W treat the two cases separately where A has semi-simple eigenvalues or a Jordan 
cell. 



4.1 When A has semi-simple zero eigenvalues 

In this subsection, we treat the case where A has semi-simple eigenvalues. Let the 
dimension of keiA be m; 

AUi = 0, (i = l,2,...,m). (4.8) 
We suppose that other eigenvalues have negative real parts; 

AUc = XaUa, (a = m + l,m + 2, ■■■,n), (4.9) 
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where ReA^ < 0. One may assume without loss of generahty that C/j's and C/q's are 
hnearly independent. 

The adjoint operator A"^ has the same eigenvalues as A has; 
A^Ui = 0, (z = l,2,...,m), 

A^ila = A;t/„, (a = m + l,m + 2,---,n). (4.10) 

Here we suppose that C/j' and UaS are linearly independent. Without loss of generality, 
one can choose the eigenvectors so that 

{Ui,U^) = = {U^,U,), (4.11) 

with 1 < i < m and m + 1 < a < n. 

We denote the projection operators by P and Q which projects onto the kernel of A and 
the space orthogonal to keiA, respectively. The projection operators can be constructed 
in terms of Ui and Ui [i = 1,2, ... , m) as follows: Let Ifp be an n x m matrix defined 

by Up = {Ui, t/2, . . . , Um) and Up by Up = (t/i, C/2, . . . , t/^), then 

P = Up{UpUp)-^Up, (4.12) 

and Q = 1 - P. 

Since we are interested in the asymptotic state as t —>■ 00, we may assume that the 
lowest-order initial value belongs to keiA: 

m 

Woito) =Y.CMU, = Wo[C]. (4.13) 

i=l 

Thus trivially, Uo{t;to) = e(*-*'')^VFo(to) = EZiCt{to)Ui. We notice that a natural 
parameterization of the invariant manifold in the lowest order Mq is given by the set of 
the integral constants C = *(Ci, C2, ■ ■ ■ , Cm) being varied. 

The first order equation ([4.5|) with the initial value VFi(to) which is not yet determined 
is formally solved to be 

u,{t;t,) = e(*-*»)^l^i(to) + f' dse^'-^^^F{u,{s;to)). (4.14) 

We remark that one may assume that the initial value VFi(to) is independent of Wo{tQ), 
namely Wiito) belongs to the Q-space, because if VFi(to) had a component belonging 
to keiA, the component could be " renormalized away" into Wq. Inserting the identity 
I = P + Q between the two functions in the integral, we have 

u,{t;to) = e^'''^>^^[W,ito) + A-'QFiWoito))] 

+ {t - to)PF{Woito)) - A-'QF{Woito)). (4.15) 
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The first term has a possibihty to give rise to a fast motion, which should be avoided 
and called "dangerous" term: The "dangerous" terms here are analogous with divergent 
terms in quantum field theory, which are subtracted away by counter terms, analogue to 
the initial values Wi here, self-consistently. Indeed it is nice that the initial value VFi(to) 
not yet determined can be chosen so as to cancel out the "dangerous" term as follows; 

W^{to) = -A-'QF{Wo{to)), (4.16) 

which satisfies PWi(to) = and is a function solely of C{to). Thus we have for the first 
order solution 

ui{t; to) = {t- to)PF - A-'QF, (4.17) 

where the argument of F is VFo[C]. Notice that Eq. (|4.17|) is consistent with (|4.16| ). 
We remark that what we have done is actually a simple thing; we have suppressed the 
unperturbed part that would be damped out as t —>■ oo. 

Now the invariant manifold is modified to Mi given by 

Mi = {u\u = Wo-eA-^QF{Wo)}. (4.18) 

If one stops to this order, the approximate solution reads 

uit; to) = Wo + e{{t- to)PF - A'^QF}. (4.19) 

Then the RG equation du/dto\tg=t = gives 

Wo{t)=ePF{Wo{t)), (4.20) 

which is reduced to an m-dimensional coupled equation, 

C,{t) = e{U„ F{Wo[C])), (z = 1, 2, ■ ■ ■ , m). (4.21) 

The global solution representing a trajectory on the invariant manifold up to this order 
is given by 

m 

u{t) = u{t;to = t) =J2ait)U^ - eA-'QFiWo[C]), (4.22) 

1=1 

with C{t) being the solution to ( |4.21|) . 

In short, we have derived the invariant manifold as the initial value represented by 
( [4.22| ) and the reduced dynamics ( |4.21| ) on it in the RG method in the first order approx- 
imation. 

The second order solution can be obtained as follows. As has been done in the first 
order case, the second order solution is formally given by 

U2{t- to) = e(*-*»'^ W2(to) + r rfse(*-^)^F'(txo(s; to))u,{s; to). (4.23) 

Jto 
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We may assume again that the initial value W2{to) belongs to the Q-space. A straight- 
forward evaluation of the integral yields 



^(t-to)A 



W2{to) - {a~'qf'a-'qf - a-^qf'pf] 



+A-^QF'A-^QF - A-'^QF'PF - (t - to) [PF'A-^QF + A-^QF'PF^ 



+ -{t-tofPF'PF, 



(4.24) 



where the argument of F and F' is VFo[C]. The initial value can be now determined so 
as to cancel out the "dangerous term", i.e., the fast moving part, as before; 



W2{to) = A~'QF'{Wo)A-'QFiWo) - A'^QF'PF, 



(4.25) 



which belongs to the Q-space. This implies that the invariant manifold is modified to M2 
represented by "W^ = VFo[C] + p[C]] p ^ eWi + e'^W2- Thus we obtain for the second 
order solution 

U2{t] to) = A-^QF'A-^QF - A-'^QF'PF - {t - to) [PF'A'^QF + A'^QF'PF] 



+ -{t-to)^PF'PF. 



(4.26) 



Notice that Eq.( 4.26| ) is consistent with ( 4.25 ). Thus the full expression of the solution 
up to the second order is given by 

u{t;to) = Wo{to) + e{{t-to)PF-A-'QF} 



1 



A-^QF'A-^QF - A-^QF'PF - (t - to) {PF'A'^QF + A'^QF'PF} 

(4.27) 



+ ^{t-to)^PF'PF 



The RG equation du/dto\ 



to=t 



reads 



Woit) - ePF - eA-^QF'Wo + [PF'A'^QF + A'^QF'PF] = 0. (4.28) 
Operating the projections P and Q to the both sides of (|4.28| ), respectively, we have 



Wo{t)-ePF + e^PF'A-^QF = 
- eA-^QF'Wo + e^A-^QF'PF = 

Firstly, we notice that the last equation ( |4.3(J| ) is reduced to 

eA-^QF'{~Wo + ePF) =0, 



0, 
0. 



(4.29) 
(4.30) 



which is identically satisfied on account of ( ^.29] ) up to this order. Thus, we end up with 
the reduced equation given by 



Wo{t) = ePF - e'^PF'A-^QF, 



(4.31) 
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which is further reduced to 



Ci = e{Ui,F-eF'A-^QF), 



(4.32) 



for i = 1, 2 m. 



The global solution giving the trajectory on the invariant manifold is given by the 
initial value as 

u{t) = W{t) = W,[C] + p[C] 

= Wo[C] - eA-^QF + e^{A-^QF'A-^QF - A'^QFTF}, (4.33) 

with C{t) being the solution to ([4.32|) . Notice that the argument of F and F' in the 
above expression is all Wq, hence the r.h.s is a function of C, i.e., u{t) = u[C]. Recall 
that C was the integral constants of the unperturbed solution. 



A couple of remarks are in order: (1) When the present formulation is applied to 
the Lorenz model |^8| around the first bifurcation point which is simple one, the result 
coincides with that given in |26]; in other words, the present formulation gives a foundation 
to the prescription adopted in that paper. (2) The present formulation using the projection 
operators and the resultant RG equation ( [4.32|) governing the slow motion resemble those 
by Mori theory [EOf for stochastic motions. 



4.2 When A has a Jordan cell 

In this subsection, we treat the case where A has a Jordan cell. We assume the Jordan 
cell is two-dimensional for simplicity, and we define the normalized vectors Ui and U2 by 

AUi = 0, AU2 = Ui. (4.34) 

The conjugate vectors Ui and U2 satisfy 

A^U2 = 0, A^ili = L/2, (4.35) 

where A"^ is the conjugate of A. The normalization condition is given by 

(C7i, C/i) = {U2, U2) = 1, (Oi, U2) = 0, (4.36) 

where ( , ) denotes the inner product. Note that {U2, Ui) = automatically holds. 

We denote by P the projection operator to the subspace (P-space) spanned by Ui 
and L/2; namely, for any vector u, 

Pu = aUi + PU2, (4.37) 

with a = {Ui, u) and 13 = {U2, u). 
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Let Q be the projection operator to the subspace (Q-space) comphment of the P-space. 
Then one can verify that 

e*^n = e*^(P + g)ti, 

= {a + (3t)Ui+ (3U2 + e'^Qu. (4.38) 

So much for the prehminaries. 

Now let us proceed to obtain the asymptotic solution to ( [4.4| ) as t — > cxd by the 
perturbation method: Since we are interested in constructing the invariant manifold, let 
us take the stationary solution as the lowest order one; 

Wo(t;to) = Co(to)C/i, (4.39) 

accordingly, the initial value reads 

Wo{h) = Co(to)C/i. (4.40) 

Notice that we have not included a component in U2 direction. The lowest order manifold 
is 

Mo = {u\u = CoUi}. (4.41) 



The first order solution is formally given by ( [4.141) . The first order initial value is 
chosen to be independent of Wq] 

Wi{to) = C,{to)U2 + QW,{to). 

A simple evaluation of the integral in ( |4.14| ) gives the first order solution as 



+{Ci{to){t - to) + apit - to) + /Split - to)^}Ui 

+{Ci{to) + Prit - to)}U2 - A-'QF + 0((t - to)""'), (4.42) 

where 

aF={UuF), (3f={U2,F). (4.43) 

The argument of F in the above expressions is W^o[Co]. The initial value can be deter- 
mined so as to cancel out the fast mode as before, namely, 

QW,{to) = -A~'QF, (4.44) 

which implies that the invariant manifold is modified to 

Ml = {u\u = CoUi + eCiC/2 - eA-^QF}. (4.45) 
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Then the solution up to the first order is obtained as 
y'{t;to) = Co{to)Ui + e 

HMt - to) + Ci(to)}C/2 - A-'QF 



{apit - to) + Ci{to){t - to) + - tofWi 



(4.46) 



The RG equation in this order is given by 

= CoU^ - e{{aF + C,)U, + {(3p - C,)U2 + jQF'CoUu 
which leads to 

Co = e{{U^,F) + C^), C^ = {U2,F). 
We now see that the trajectory in the invariant manifold Mi is given by 



(4.47) 



(4.48) 



n(t) 



W{t)c^Wo[Co] + Wi[Co], 
Co(t)C7i + eCi(t)C/2 - eA-^QF, 



(4.49) 



with Co(t) and Ci(t) being governed by ( [4.48| ). Notice that n(t) is a functional of Co(t) 
and Ci(t). 

The second order solution is obtained, similarly. The formal solution to the second 
order equation is given by Eq.( |4.23| ). Let g = Ci{to)Ui + apUi + PfU2 and h = 

——QF + CiU2- Then, a simple manipulation as before gives 

Mt;to) = e^'-"^^^[W2ito) + {A-'QF'h + A-^QF'g}] 

+ (t - to) {-A-^QF'g + PF'h} - {A-^QF'h + A~^QF'g} 
+0((t-to)2). (4.50) 



Thus the initial value is determined to be 

W2(to) = - (A-^QF'h + A-^QF'g 



(4.51) 

so that the fast modes disappear; notice that W2 belongs to the Q-space. The invariant 
manifold is modified now in an apparent way; so we do not write the specification of the 
second order manifold M2. Hence we have 



U2{t; to) = (t - to) [-A-^QF'g + PF'h] - {A-^QF'h + A-^QF'g} 
+0{{t-t,f). 



(4.52) 



Applying the RG equation to n = + ei^i + £^'"■2 thus obtained, we have 

CoC/i - e{aF + Ci)Ui + e{Ci - (3f)U2 - eA~^QF'CoUi 
~e^{-A-^QF'g + PF'h}-e^[A~^QF'CiU2 + A~^QF'CiUi} = (4.53) 
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Operating P and Q to ( [4. 531 ), we obtain 



= CoUi + eCiU2-eg -e'^PF'h, (4.54) 
= ~e^A-^QF'CiUi-eA-^QF' {CoUi + eCiU2-eg}. (4.55) 

We remark that from ( ^.54| ) and ( [4.55| ), 

QF'Ui = (4.56) 

must hold as a compatibihty condition, which gives a necessary condition for the RG 
method to work for higher approximations. With this condition taken for granted, g is 
reduced to = [3pU2- 

Equating the components in the Ui, U2 in ( [4.54| ), we have the reduced dynamics as 
follows, 

Co = e{{lJ^,F + eF'h) + C^), (4.57) 
Ci = {U2,F + eF'h), (4.58) 

where the argument of F and F' is Wo{t) = Co{t)Ui. 

The trajectory on the manifold M2 is given by 

u{t) = Wit) = Woit) + eWiit)+e''W2it), 
= Co{t)Ui + eCi{t)U2-eA-^QF 

-e^A~'QF'h + {U2, F)A~^QF'U2}. (4.59) 

Comments are in order: The solution is solely described by the coordinates Cq and Ci 
representing the P space as in the previous subsection. Conversely, the dynamics can not 
be described only by the coordinate in the zero-th manifold in this case; the dimension of 
the invariant manifold is increased from that of the unperturbed invariant manifold. We 
remark that the formulae obtained above are completely consistent with those given for 
the Takens equation in §3. 



4.3 Practical way 



We have formulated the RG method so that it is clarified that the method concerns with 
the initial values and the invariant manifold is constructed as the initial value perturba- 
tively: The initial values are determined so that terms proportional to the unperturbed 
solution and those representing fast motions disappear in the perturbed solutions. In 
effect, the special solutions of the higher order equations are composed of secular terms 
which are proportional to the unperturbed solution and vanish at t = to and solutions 
independent of the unperturbed solution. We remark that this way of construction of the 
perturbative special solutions have been adopted in [^, ^]. In this subsection, having 
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known the above fact, we present the rules for constructing the special solutions in the 
form of an operator method; a detailed account of this method is given in Appendix A. 
This subsection will constitute a practical summary of the results obtained in the previous 
subsections. 

When we try to obtain an invariant manifold and the reduced dynamics on it, the 
solution to the unperturbed equation ([4.4|) is given by a stationary one 



no(t;to) = VFo(to). (4.60) 
Then the first order solution is given by 

niit;to) = ^^F(no) = ^^(PFK) + QFM), 

= {t~to)PF+^QF. (4.61) 



Here ( |A.16| ) and (|A.17|) have been used and A ^ is written as ^ 



Similarly, the second order solution is given by 



^ iP + Q)F'{uo){{t-to)PF+^QF}, 



dt-A' ' ' "'^^ -A 

i(t - ufPF'PF + (t - t^)PF'^QF 

-{{t + ^}QF'PF + J^QF'^QF, (4.62) 



which coincides with ([4.26|) . Here we have used the formulae (|A.20|) , (|A.21| ) and (|A.22| ). 
The efficiency of the operator method is apparent. 

Next, let us consider the case where A has a two-dimensional Jordan cell. We take a 
stationary solution as the zeroth-order one; 

no(t;to) = Co(to)C/i. (4.63) 



Notice that the kernel of A is not yet fully spanned by the solution. Then the first 
order solution is given by a sum of the remaining component of kerA and the special 
solution; 

= Ci(to)(t - to)C/i + Ci(to)C/2 + -^-jiPFiuo) + QF{uo)), 
= Cr{t,){t-t^)U^ + Cr{to)U2 

+{it - to)PF +^{t- to f{if2, F)U,} + -^QF, (4.64) 
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which is found to coincide with Eq.( |4.42| ). Here ( |A.17| ) and ( |A.27| ) have been used. The 



unperturbed solution ^^^^{O} may be obtained in the following way: Putting ^^^^{O} 
a{t)Ui + b(t)U2, one has a coupled equation d = b, 6 = 0, which is solved to yield 
a(t)=Ci(to)(t-to), bit) = CM. 

The second order solution is given by 

U2{t;to) = g ]_ ^ F'{uo)ui{t;to), 

^ {P + Q)F'{uo)u^{t;to). (4.65) 



dt-A 



Then the r.h.s. is a sum of terms which are in the form — tQ)'^PG{uQ) or ^^^^ (t — 

to)"'Q(?(wo)) which are calculated in Appendix A. Actually, since terms proportional to 
(t — to)" ^ 2) do not contribute to the RG equation nor to the resulting trajectory, 
one needs not calculate the terms of the form ^^^(t — to)"'PG{uo) with n > 1. Thus we 
easily reach the final result given in ([4.52|) . 



4.4 Normal form 



As a simple example, we try to reduce the following two-dimensional evolution equation 
with a Jordan cell; 

X = y + aix^ + bixy + Ci?/^, ij = a2X^ + b2xy + C2?/^, (4.66) 

with Oj, bi and Cj {i = 1,2) being constant. As in the Takens equation examined in §3, 
we make a scale transformation 

x = eX, y = eV. (4.67) 

Then the equation is reduced to 

{dt- A)u = eF{u), u= *(X,F), (4.68) 

where 

Uoj' ^^''^ [a2X^ + b2XY + C2Y^)- ^^-^^^ 
Here the P-space is spanned by 

U,= Q, and U2= (^^, (4.70) 
which satisfy AUi = and AU2 = Ui. And tJi = Ui {i = 1,2). 
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Expanding as u = Uq + eui + ■ ■ ■ , with Ui = Yi), we first solve the equation 

around t ^ to with the initial value W(to) = Wo(to) + eWi(to) + ■ ■ ■■ 

We take the stationary solution as the lowest order 

Uo{t;to) = Aoito)Ui, (4.71) 

where v4o(to) is an integral constant. Accordingly, VFo(to) = ^(^ol^o), 0); the unperturbed 
manifold is the X axis. 



According to the general argument given in §4.2, we need to calculate the following 
quantities; 

ap = {Uu Fiuo)) = aiAlito), Pf = {U2, Fiuo)) = a2Al{to). (4.72) 
Notice that QF{uq) = 0, identically. Then one has 

ni(t;to) = {Ai(to)(l + aiAi(to))(t-to) + Y^o(^o)(t-to)'}C/i 

+{Mto) + a2Al{to){t - to)}t/2, (4.73) 

where Ai{to) is another integral constant. Accordingly, Wi{to) = *(0, Ai(to))- Then 

M, = {u\u= \Ao,eA,)}. (4.74) 

If we stop at this order, the RG equation reads 

Ao = e{Ai + aiAl), ^1 = 02^2. (4.75) 

This is the Takens equation; notice the trajectory is given by x(t) = eX{t] t) = eAo{t), y(t) - 
eY{t; t) = e^Ai{t). It means that the RG method gives mechanically the normal form |Q 
of the reduced equation. This is confirmed by proceeding to the second order. 

In the present case, g = and h = Ai{t 0)1/2, and 

PF'{uo)h = Ai{biAoUi + (4.76) 
Then the second order solution is found to be 

U2{t; to) = {t- to)Ai{hAoUi + 62^01/2} + 0{{t - to)""'), (4.77) 
Accordingly, W2(to) = 0. 

Up to this order, we have u = Uo + eui + e'^U2. Then the RG equation reads 

io = eAi + e{aiAl + ehAoA^), 

Ai = 02^^ + 662^0^1, (4.78) 



which is the normal form when the linear matrix is of Jordan form|Q. The trajectory 
is given by x{t) = eAo(t), y{t) = e^Ai(t). Thus one sees that the RG equation gives the 
normal form of the reduced evolution equation on the invariant manifold. We remark 



that when ai = bi = 0, the RG equation is nothing but the Bogdanov equation|51 . 



We remark that one needs the second order solution to give the Bogdanov equation, 
while the Takens equation is obtained in the first approximation. 
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4.5 An extended Takens equation 



As the final example, we deal with an extension of the Takens equation to a system with 
three-degrees of freedom; 



X 



y + eax , y 



ebx 



z = -z + ef{x,y,z), 



(4.79) 



where f{x,y,z) is analytic function of {x,y,z). We shall show the compatibility condi- 
tion ([4.56|) , which becomes relevant only when the system has more than two-degrees of 
freedom, gives a restriction to the form of f{x,y, z). 



With u = ^{x, y, z), ( |4.79| ) is converted to 



where 







A)u = eF{u), 




/ 1 \ 




1 ax"^ \ 


A = 





, F{u) = 


bx"^ 




loo -ij 




\fix,y,z)J 



Notice that A has a two-dimensional Jordan cell; 



AU, = 0, AU2 = C7i, AU3 



(4.80) 



(4.81) 



(4.82) 



where Ui = *(l,0,0),C/2 = *(0, l,0),C/3 = *(0,0,1). The projection operator to the 
subspace {Ui,U2} is given by P = diag(l,l,0), while Q = 1 — P = diag(0,0, 1). We 
also notice that Ui = Ui {i = 1,2) in this simple example. We are interested in the 
asymptotic behavior of the solution at t — >■ 00. We first solve ( [4.80|) around t to with 
the initial value W^(to) at t = to by the perturbation theory. The solution may be written 
as u{t] to, W {to)) , which is expanded as tt = uq + eui + e'^U2 + ■ ■ ■ . The initial value which 
is determined self-consistently with u is also expanded a.s W = Wq + eWi + e^W2 + ■ ■ ■ . 



When t — > 00, we may take the stationary solution as an asymptotic one to the lowest 
order equation; 

Wo(t;to) = Co(to)C/i. (4.83) 

Accordingly the initial value reads VFo(to) = Co(^o)C^i, which implies that the unper- 
turbed invariant manifold is given by 

Mo = {u\u = CoC/i}, (4.84) 

namely the x axis. 

According to the general formulation given in the previous sections, to obtain the first 
order solution, we only have to evaluate 

ap = {U,,F{uo))=aC^,{to), Pf = {U2, F{uo)) = bC'^ito), 
A-'QF{uo) = -/(no)t/3. (4.85) 
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Then we end up with 

ui{t; to) = {aC^itQ)it - to) + Ci(io)(i - ^o) + lbC^ito){t - to)''}Ui 

+{bC'o{to){t - to) + C\{to)}U2 + f{uo)U^. (4.86) 

Accordingly, 

W^{to) = C^{to)U2 + /(wo)C/3, (4.87) 
which imphes that the modified invariant manifold is given by 

Ml = {u\u = \Co, Ci, /(Co, 0, 0))}. (4.88) 

If we stop at this order, the full solution is given by w ~ Uq + ewi- Applying the RG 
equation, we have 

Co = e{aCl + Ci), Ci = hCl (4.89) 
And the trajectory on the manifold M~Mi is given by 

u{t) = Co{t)Ui + eCi{t)U2 + ef{Co{t), 0, 0)Us. (4.90) 

To go to the second order, we first need to examine the compatibility condition; 

df{x, y, z) 



= QF\uo)Ui = 



dx 



Us, (4.91) 

U=Uo 



with Uo — *(Co, 0, 0), which means that when y — z — 0, f{x, y, z) does not depend on 

X. 

With this condition assumed, we can proceed to the second order. To obtain the second 
order RG equation, we notice the following; h — —A~^QF{uo) + CiU2 = f{uo)U3 + CiU2 
and hence F'{uo)h — 0. Thus the RG equation which gives the evolution equation of 
Co,i(i) is not modified by the second order perturbation. 

To obtain the second order correction to the trajectory, we need to evaluate the fol- 
lowing; 

-A-'QF'h = (f(^o)% + Ci^)Us: 

-A-'QF'U2 = -|^^3, (4.92) 

where the derivatives are evaluated at u = tto = *(Co(t), 0, 0). Thus the second order 
correction of the initial value reads 

W2{to) = ifiuo)^ + C^^ - Pf^)Us. (4.93) 
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Hence the trajectory in the second order approximation is given by 



u{t) = W(t) = C„{t)Ui+€Ci{t)U2 + ef(ua)U, 

+^^(/M| + C.,|-,3.|)t/3. (4.94) 

Here Co,i(t) are governed by the Takens equation ( |4.89| ). We see that the higher order 
terms does not affect the dynamics but modifies the trajectory only in the U2 and the 
Q-direction. 

Similar discussions can be made for the following extended Bogdanov equation, 

X = y, y = ax'^ + bxy, z = —z + ef{x,y, z), (4.95) 
where f{x,y,z) is an analytic function. 



4.6 Discussion 



We have applied the naive perturbative expansion as the starting point of the RG method. 
However, the naive perturbation expansion is not always a good starting point. There are 
cases where a scaling transformation is needed to convert the unperturbed equation to a 
non- Jordan form before applying the perturbative expansion. Let us take the following 
example; 

^=(^}~n')^^ (4-96) 



e . 

with u = ^{x,y). The exact solution reads 

u{t) = Ae^'U+ + Be-^'U-, (4.97) 



where A = y'(l - e)e and U± = *(1, ±^e/(l - e)); i.e., 

x{t) = Ae^' + Be-^\ y{t) = ^e/{l - e) ■ (Ae^* - Ee"^*). (4.98) 



Due to the appearance of the singular term y/e, the naive perturbation does not give 
a sensible result even in the RG method. In this case, we first try to convert the equation 
so that the converted equation has no Jordan nature in the unperturbed part. This can 
be performed by a scale transformation; x = e"^, y = e^r], t = e^r. To make the 
unperturbed part to be a non- Jordan form, we choose that [3 — a + u = l + a — /5 + z/ = 0, 
which is satisfied with a = 0, [3 = —v = 1/2, i.e., x = ^, y = y/er], t = r/y/e. The 
converted equation reads 

f g^e, (4.99) 
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which can be now solved by the RG method with no difficulty. The result obtained by 
the RG method up to O(e^) reads 

x{t) = Ae^'* + Be-^'\ y{t) = ^{1 + e/2)(Ae^'* - Be"^'*), (4.100) 

with A' = ^(1 -e/2). 
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5 Applications I 



In this section, we present examples of non-linear equations for which the unperturbed 
linear operator has no zero eigenvalues but a pair of eigenvalues Aj (i = 1,2); (i) Ai,2 = 



and (ii) Ai,2 = ±A where uj and A are real numbers. The first case is discussed in [g5|, p6 
as well as in §3. We present it here for completeness. 



5.1 Brusselator 



This is an example of systems showing a Hopf bifurcation. An RG treatment of generic 
systems with a bifurcation has been given in where an emphasis is put on the relation 
of the RG method with the envelope theory. Here we treat this interesting example in 
the present formulation emphasizing the aspect of the RG method as the one to construct 
attractive manifolds. 



The Brusselator is given by 

- = BX-X^y + Dy^, (5.1) 

where A{> 0),B{> 0),Dx and Dy are constant. We here treat a uniform system, hence 
the terms with the spatial derivatives vanish. The steady state is given by (Xq, Yq) = 
{A, B/A). Shifting the variables as ^ = X — Xq, r] = Y — Yq, and defining u = *(,^, r/), 
we have 

with 

f{^,r])=B/A-e + 2A^r^ + eV- (5.3) 

The linear stability analysis shows that when B exceeds the critical value B^ = 1 + A"^, 
there arises a bifurcation. 



Now let us analyze the slow motion and the slow manifold around the bifurcation 
(critical) point. We define the following variables 

/i = (5 - B^)/Bc, e = •y/j/i|,and x = sgn(/i), (5.4) 

accordingly, fi = xe^. We first expand u and the initial value u{t = toS^o) = W(to) as 
Taylor series w.r.t e: u = eui + e^U2 + e^ns + ■ ■ ■, and W = eWi + e^W2 + e'^Vl^s + ■ ■ ■. 
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The first order equation reads 



with 



{dt - Lo)wi = 0, 



A" A" 



-(^2 + 1) - Ay- 
The solution is readily found to be 

u{t; to) ^ C{to)Ue''^^ + c.c, 

with uj = A and 



U 



1 

■ 1+-L4 



Here C is the (complex) integral constant. Accordingly, 

Wi{to)^C{to)Ue"^''' +C.C.. 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



A simple manipulation using the formulae given in Appendix A gives the higher order 
terms as follows; 



U2{t;to) = {CV+e'""* + c.c.} + |C|Vo, 



where 



V 



+ 



^2-1/0 



3A3 \l + 2iAj 



A3 V 1 y ' 



(5.10) 



(5.11) 



and 



where 



^{{t- to)U + ^U*}e'^' + ^{U + )-U*)& 

IlUJ 4:IUJ 2 



2 



* \ _3ia;t 



+C.C., (5.12) 



Ci = xBcC+{ T^ + i 



A^ 3^3 
Cs = {2B,V+^/A + 2A{V+r, + V+^U, + U,)}C^ 



}\c\'c, 



(5.13) 
(5.14) 



with V^^ being the ^-component of V+ and so on. Here the initial values have been 
chosen to be 



W2ito) = {C^V+e^'^'^ + c.c.} + \CfVo, 
Ci ... C, 1. 



= 

Here C7* is the complex conjugate of U. 



Aicu Aicu^ 2 ' 



+C.C.. 



(5.15) 
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Collecting all the terms thus obtained to have the approximate u{t; to) and applying 
the RG equation to it, we have dC/dt — e^Ci = 0, or 



with 



2 + ^2 -AA^ + 7^2-4 
+ 1 — 



The attractive manifold is given by the initial value as 

u(t) = W{t) ^eWi{t) + e^W2{t) + e^W{t), 



e{C(t)C/e^"* + c.c.} + e'[{C'{t)V+e''^' + c.c.) + \C\'Vo] 

I ^3 



+ c.c. 



(5.16) 



(5.17) 



(5.1^ 



These results coincide with those obtained in the reductive perturbation method[18 



5.2 Unstable motion in Lotka— Volterra system 



The Lotka- Volterra (LV) equation |5^ is known to be integrable, although the exact ana- 
lytic solutions of it are not known. The equation was already treated in the RG method 
by one of the present authors (TK), and an approximate solution was constructed explic- 
itly around the non-trivial fixed point A numerical comparison of the results with 



the exact solution shows that the solution given by the RG method well approximates 
the exact solution in a global domain even when the small parameter in the equation e, e' 
(see below) are as large as 0.8 |Q. The purpose of the present subsection is to apply the 



RG method for analyzing the equation around the unstable fixed point (the origin). 
Lotka- Volterra equation reads 

= ax-exy 
- —by -\r e xy 

Here, a; = x{t),y = y{t) and a, 6, e, e' are positive constants. In this work, we treat the case 
where < e < 1 and < e' < 1 so that the perturbation theory can be applied. 

There are two fixed points; (i) x = ?/ = and (ii) x = h/e' , y = a/e. One can see 
that the fixed point (i) is asymptotically unstable. An approximate but globally valid 
solution around the second fixed point was obtained in [^. In this paper, we treat the 
first fixed point. 

With the new variables defined by x = ^S,, y = V ( p.l9| ) is converted to 

^ = Au + eF{u), (5.20) 
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where 



"=0. •4=(o -i). F(u) = -Ol[\). (5.21) 

One see that the fixed point (i) is asymptotically unstable in the linear approximation. 
Our aim here is to obtain the dynamics around the unstable fixed point and construct 
the unstable manifold using the RG method. 

For later convenience, we introduce the normalized eigenvectors of A; 

AUi = aUi, AU2 = -bU2, (5.22) 
where Ui = *(1,0) and U2 = *(0,1). 

Now let us apply the RG method to construct an approximate solution valid in a global 
domain. We suppose that the equation is solved up to arbitrary t = to from the time, 
say if: = 0, at which the genuine initial condition is imposed. With this up-to-dated initial 
value W{to), we try to construct the solution to ( |5.20|) around t ~ to by the perturbation 
theory, expanding u a.s u = Uo+eui + e'^U2 + 0{e'^). Thus u may be written asu = u{t] to). 
We also expand the initial value W(to) = Wq + eWi + t^W2 + o{e'). 



The equations to be solved are 

(^-A)'ao = 0, (5.23) 

(^-A)ni = -^or7o((7i-(72), (5.24) 

{j^-A)u2 = -{^om+^iVo){Ui-U2), (5.25) 

and so on, where = Ui {i = 0, 1,2...). It turns out that one needs to treat sepa- 

rately depending on whether a ^ h 01 a = h. 



Case A: a = b 

This is the interesting case where secular terms appear and the RG method plays a role 
to construct an approximate solution in a global domain. The unperturbed solution reads 

uo = Ci{to)e^'Ui + C2(to)e~"*C/2, (5.26) 

implying the initial condition VFo(to) = Ci{to)e"-^°Ui + C2(to)e~"*°C/2. The first order 
solution now reads 

Ui = -C,C2jr^iUi-U2) = ^iUi + U2), 

Ot — A a 
= W.ito), (5.27) 

which is constant and independent of Uo{t). 
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Using the formulae in Appendix A, the second order solution is also obtained. Col- 
lecting all the terms thus obtained, we have an approximate solution valid around t ~ to; 

u{t:to) = Ci(to)e'^*C/i + C2(to)e-'^*C/2 + e^(C/i + i72) 

a 



a 2a^ 
+ {^e»* + ^(t-to)e-'^*}C/2]. (5.28) 



Applying the RG equation du/dto\tQ=t = to (|5.28|) , we have 



C<^ + e2^ = 0, ^2-6^^ = 0. (5.29) 



a a 
Noting the constraint C1C2 = const. = c given by (|5.29|) , we reach 



Ciit) = C2(t) = -e^\ (5.30) 

Cl 

with Cl being a constant. 

Our solution is given by the initial value W{t)] 

u{t;t) = W{t) = Ci(t)e'^*C/i + C2(t)e^'^*C/2 + e-(C/i + C/2) 

a 

with Ci^2if) given by ( |5.30| ). Introducing C2 by c = C1C2, the respective components are 
given by 

mt) = Cie('^-AP)* + !^ + !!Me-(-^)* + o(e^), (5.32) 

Vit;t) = C2e-(^-^)* + ^ + 4^e(^-^)* + o(6^). (5.33) 
Here, the constants Ci and C2 are determined by the initial condition imposed at t = 0. 



Some remarks are in order: (1) Our solution ( |5.31| ) shows that the speed to approach 
to and to escape from the origin is shifted from a to a — e^ciC2/a = a, which is dependent 
on the initial condition as expressed by ci and C2. (2) The unstable manifold can be 
constructed explicitly from (|5.32|) and (|5.33| ). Solving e"* and e~"* from these equations 
and making the product of them, we have 

c= {X -eY){Y -ex), (5.34) 



where X = — xo, Y = y — yo, & = and xo = yo = • ( |5.34] ) shows that the 



unstable manifold is the hyperbolic curve which has the non-orthogonal asymptotic lines 
coming out of the point *(xo,a;o) with the slopes 9 and 1/9, respectively . 
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Case B: a ^ b 

This is a rather trivial case because no secular terms appear, and the RG method does 
not play any role. Therefore we only write down the result: 

+ 

rj{t) = 
+ 

where Ci and C2 are constant. 



262 

2 



(5.35) 



6*26 + e 



2g(a-6)t 



a 



za"' o(a — 0) 



(5.36) 
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6 Applications II: Pulse interactions 



In this section, we apply the RG method to obtain the dynamics of interacting pulses 
(or fronts) in one dimension. Some yeas ago, a systematic method was developed by Ei 
and Ohta[45| to this problem on the basis of the phase dynamics approach which 



involves a solvability condition. We shall show how the dynamics of interacting pulses 
obtained by them can be derived mechanically in the present method virtually without 
any assumptions. In other words, we shall derive the phase equations describing the front 
and pulse interactions in the RG method for the first time, although interface dynamics in 
spinodal decomposition [2^, a diffusion equation|30] and Kuramoto-Sivashinski equation 



as the phase equations have been derived in the RG method by others. As repre- 
sentative examples of pulse dynamics, we take up the kink-anti-kink interactions in the 
time-dependent Ginzburg-Landau (TDGL) equation and the soliton-soliton interaction in 
the Kortweg-de Vries (KdV) equation. 



6.1 Kink-anti-kink interaction in TDGL equation 



The TDGL equation we study is given by 



(6.1) 



with 



m = lu{l-u'), (6.2) 

where u = u{x,t) is a real scalar function and < e < 1. We remark that (|6.1| ) has a 
stationary solution 



X — h 

u{x, t) = ± tanh = ±U{x — h), (6-3) 

where U{x) satisfies U{±oo) = ±1 and f/(0) =0. U{x — h) {—U{x — h)) is called a kink 
(an anti-kink) at the position x = h. 



6.1.1 Interaction of one kink and anti-kink 

We first consider the interaction of one kink and anti-kink. Suppose the initial data at 
t = are given by 

Uoix; xi(0), X2(0)) = U{x - xi(0)) + {-U{x - X2(0)) - 1}, 
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where X2(0) — Xi(0) >> e; i.e., the kink (anti-kink) is located at x = Xi(0) (x = X2(0)). 
Since Uo{x; xi(0), X2(0)) is not a solution to Eq. (|6.1|) , the position of the kink xi and the 
anti-kink X2 will move slowly. Our task is to find the equation governing the dynamics of 

Xi{t) and X2{t) at t > 0. 

To apply the RG method to solve this problem, we first identify the small parameter 
and the integral constants in the unperturbed solution, which will move by the perturba- 
tion. When X2 — Xi » e, Uo{x; Xi, X2) in the neighborhood of x = Xi may be represented 
by 

Uo{x; Xi, X2) = Uo(t; Xi){x) + 6s{x — X2), Uo(t; Xi){x) = U{x — Xi), 

where (5 is a small parameter because the effect of the anti-kink —U{x — X2) — 1 is small 
around x = Xi; the order of 6 is exp( — (x2 — Xi)/e). This implies that as the unperturbed 
solution around x ~ Xi, we may take uo{t;xi), where Xi is a constant. Thus we are led 
to represent the solution at t ~ Vto by 

u{x, t) = uo{t; xi{to)){x) + 6s{x - X2{to)) + v{x, t), (6.4) 

with v{x,t) as small as 6. Inserting ( |6.4| ) into (|6.1|) , we have 

dtv = Av + if'iuo) - ni))6s + 0{6^ + \v\^), (6.5) 

with 

Av^(e^^ + f'iuo)y (6.6) 

Here, we have made use of the fact that uq and 6s + 1 are stationary solutions to Eq. (|6.1|) ; 
we also note that /'(I) = /'(—I). 



(|0| ) is an equation of the type discussed in §§4.1. Indeed, owing to the translational 
invariance of the TDGL equation ( |6.1| ), the self-adjoint operator A has a zero eigenvalue 
with the eigenfunction Ui together with the corresponding adjoint eigenfunction Ui 

Ui{x)=U,{x)=dMix-x,). (6.7) 

Therefore, one can obtain the approximate solutions of ( |6.5| ) and the dynamics of xi{t) 
according to the procedure developed in §§4.1. 

The solution of the eigenvalue problem for the operator A may be seen in text books 
on quantum mechanics|^^. Let P be the projection operator onto the kernel of A and Q 



the one onto the subspace compliment to the kernel A; we call the respective subspaces 
P- and Q-space. We notice that 

Pu=^^I^U,ix), (6.8) 

< Ui,Ui > 

where < u,v > denotes the inner product defined by 

<u^v>= / u{x)v{x)dx. (6.9) 



45 



Now let us apply the RG method to solve f, thereby obtain the dynamics of xi{t). 
First we expand v as 

v = dvi + 0{6^). (6.10) 

The equation for vi read 

dtv, = Av^ + if'iuo) - 

In this case, VFo(^o) in §§4.1 is U {x—xi{to)). Wi(to) in Q-space is given by {—A)^^Q{f'{uo) — 
f'{l)}s by (|4l6|) and vi is 



v,{t; to){x) = (t - to)P(/'K) - + ^-'Q(/'M - f'{l))s. (6.11) 

Thus, we have an approximate function 

u{t; to){x) = U{x - Xi(to)) + Svi{t; to){x) + 5s{x - X2(to))- 
Applying the RG equation to u{t; tQ){x), one has 

= -XAU{X - Xi) - P{f{uo) - f'{ms{x - X2) + 0(5^), 







du 



to=t 



leading to the dynamics governing xi, 

<U^,if'iuo)-nl))6s> 



Xi 



+ 0(5^ 



<Ui,Ui> 

< d^jx - Xi), {f'iUjx - Xi)) - f'{l)){-U{x - X2) - 1} > 

<u,,u,> 

< d,U{x), U\U{x)) - /'(l)){-f/(x - (X2 - Xi)) - 1} > 



+ 0(5^ 



<f/i,t/i> 



12ee" 



+ 0(5^ 



(6.12) 



which coincides with the result by Carr-PegoQ], Fuco-Hale[0 and Ei-Ohta^^ where 
an explicit evaluation of the inner products is also given. 



6.1.2 Kink-anti-kink interaction in the presence of infinite kinks and anti- 
kinks 



We next consider an initial value problem of Eq.( |6.1|) with the initial condition where 
infinite kinks (anti- kinks) are located periodically a,t x = h + 2n {x = —h + 2n) with 
n = 0, ±1,±2... PI , We suppose that the intervals of a kink and the neighboring 
anti-kinks are much larger than the width of a kink (anti-kink), i.e., h ^ e. In this 
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situation, the problem may be formulated as an initial value problem with one kink at 
X = h in a finite domain < a: < 1 with a Neumann boundary condition, i.e.. 



du 
dx 



du 

x=o dx 



(6.13) 



x=l 



The initial profile in this case can be approximately represented by the function 

Uo{x;h) = U{x-h) + {-U{x + h-2)-l} + {-U{x + h) + l}, (6.14) 
= uo{x)+5ir{x + h-2) + 52i{x + h) (6.15) 

where the second and the third terms denote the small effects coming from the anti-kink 
at X = 2 — h and x = —h, respectively; the smallness of the effects are represented by 
the parameters 6i and 62, which are the orders of exp(— 2/i/e) and exp(— 2(1 — h)/e), 
respectively. 

Let us represent the solution for t > by 

u{x,t) = Uo{x] h) + v{x,t), (6.16) 

where [^(x, t)| is supposed to be as small as 61 and 62. Inserting ( p.l6| ) into Eq.(|6.1D, we 
have 

dtv = Av + if'iuo) - /'(l))(5ir + 62i) + 0{Sl + 5l + \v\^) (6.17) 

in a similar manner with (|6.5D. Here, Av = (^^^ + f'i'^o))^ and the eigenfunction 
associated with zero eigenvalue in this case is Ui = dxU{x — h). 

Let us apply the RG method to solve v, thereby obtain the dynamics of h{t). First 
we expand v as 

V = Sivi^o + <^2%i + 0{6l + Si). (6.18) 

The equations for t>i and fo,i read 

dtVi,o = Avi,o + (/'(«o) - /'(l))r, (6.19) 
dtVo,i = Avo,i + ifiuo) - (6.20) 

respectively. These are of the same form as treated in §§4.1, hence readily solved to be 
Vi,oit; to)ix) = {t- t,)P{f{u,) - /'(l))r + A-^Q{f{u,) - /'(l))r, (6.21) 

where r = r{x + h — 2) and h = h{to), by choosing the initial value 

Vi,o{to; to){x) = -A-'Q{f'{uo) - ni)Hx + h{to) - 2), (6.22) 

and fo,i(x,t) with r i. Thus, we have 

u{t;to){x) = U{x- h{to)) + 5iVifl{t;to){x) + 52Vo^i{t;tQ){x) + Sir + S2L 
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dto 



to=t 



Applying the RG equation to u{t; to) (a;), one has 

-hU'ix -h)~ Pif'iuo) - /'(l))((5ir(x + h-2) + 52(i{x + h)) 
+0{5l + 5l), (6.23) 
leading to the dynamics governing 

<[/i,(/'M-f(l))(5ir + 52^)> 



— e- 



2h 2{l-h) 1 „ , n o 

-12e ( e"- - e — \ + 0{5l + 5? 



2>i 



which also coincides with the known result by Carr-Pego[44] and Fusco-Hale 46 . 



(6.24) 



6.2 Soliton-soliton interaction in KdV equation 

The KdV equation reads 

dtu + Qud^u + dlu = 0, 
which has a one-pulse solution given by 

' \fc{x — ct) 



u{x, t) = -sech^ 
2 



= (p{x — ct; c), 



with c being a velocity. 



(6.25) 



(6.26) 



We consider the following problem; When two sufficiently separated pulses with almost 
the same velocities are located at Xi(0) and ^2(0) where 1x2(0) — Xi(0)| S> 1/a/c, how will 



the locations Xi{t) {i = 1,2) change at t > 0? To solve this problem, Ei and Ohta 
started with the ansatz 



u{x, t) = ip{x — ct — Xi; c + Xi) + ip{x — ct — X2; c + X2) + b{x — ct, t). 



(6.27) 



In the present work, we shall apply the RG method without any ansatz to this problem 
and show that the same evolution equation as that obtained by Ei and Ohta is derived, 
thereby give a foundation of their treatment. 

To study the problem, it is convenient to change the independent variables to 

t = t, z = x-ct, (6.28) 

namely to change to the co-moving frame with the pulse. Then the equation is converted 
to 



dtu + F[u] = 0, 



(6.29) 
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with 

F[u] = -cd,u + Qud.u + d^u. (6.30) 

We remark that 

F[ip{z-b;c)]=0, (6.31) 
with b being an arbitrary constant. 

Let us suppose that the solution around t = to > is given by 

u{z,t) = (p{z - zi{to)]c) + ip{z - Z2{to);c) + v{z,t), (6.32) 
= (^(1) + (^(2) + (6.33) 

where \z2 — Zi\ is sufficiently large. To study the effect coming from the other pulse, it is 
sufficient to consider the case of either z ^ zi ot z Z2. Then v is considered to be small. 



Substituting Eq.( |6.32|) into (|6.29| ), we have the equation governing v as 



dtv + + v?(2)]i; + 69,(<^(iV^'^) + 0{\v\^) = 0. (6.34) 

Here we have used the identity F[ip^^'^ + ^(2)] = f\^^W + (^(2)] _ - ^ 

69^((^(iV^^^) on account of Eq.(|63TD. 



Now let z ^ Zi, then (^^"^^ is small and we may put 

(^(2) = - Z2(to)) = 5g{z - Z2{to)), (6.35) 
with a small parameter 6 whose order is e~^^^^^^'^\ Then, (|6.34|) of v becomes 

dtv = - 66d,{ip'^^^g) + 0{6^ + \v\^), (6.36) 

where 

Transforming z — -2i(to) to 2;', we see ( |6.36| ) becomes 

dtV = Av- 5Qd,,{^g^^^) + 0(5^ + \v\^), (6.37) 
where A = -F'{ip) and 6g'^''\z') = ip{z' - h) with h = Z2(to) - ziito). 



(16.371 ) is an equation of the type discussed in §§4.3; the hnear operator A has a Jordan 



cell reflecting the translational invariance of KdV equation and the arbitrariness of the 
velocity of a pulse, 

AUi = 0, AU2 = Ui, (6.38) 

where 

Ui = d,^, U2 = -d,^. (6.39) 
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The adjoint operator A'' of A reads 

At = -cd,> + dl + 6(^a,,, (6.40) 
which has also a two-dimensional Jordan cell. The zero mode U2 of A'^ is found to be |^ 



U2{z)=v{z';c) (6.41) 



while these eigenf unctions have not been normalized as ( [4.36| ) yet. It is known |^8| that 
there exists a function Ui which satisfies 

A^Ui = U2. (6.42) 

The explicit form of it will be given later. 

Let P be the projection operator onto the subspace spanned by f/i and U2, and Q 
onto the subspace compliment of the P-space. 

Noting that \v\ < 0{6) and expanding v as 

v = 5vi + 5^2 + ■■■, (6.43) 

one has for vi from (|6.37| ) 

dtVi = Av,-6d,,{^g^''^). (6.44) 

Now we can proceed according to the general procedure given in §§4.3 as follows: Note 
that Uq and hence l^o('^o) in §§4.3 correspond to ip{z — zi(to);c) because both of them 
comes from the translation invariance with respect to z. Let us solve the equation ( p.44| ) 
with the initial value Wi{to). Then, we have 

v^{t;to) = e(*-*°)^W^i(to) - 6 [\^'-'^^d,iipg^^^)ds, (6.45) 

with 

WM = C{to)U2 + QW^{to) (6.46) 

as discussed in §§4.3. Here we note that although 1^1(^0) could have a Ui component, its 
effect can be taken into account by a redefinition of zi{to) which is not yet determined. 
According to the discussion given in §§4.3, the Q-component of Wi should be 

QWiito) = -A~'Q{-6dAvg^''^))- (6.47) 

Thus we have 

Vi{t;to){z') = !^{t-to)C{to) + {t~to)a + ^{t-toY(3^Ui 

+{C(to) + {t- to)P}U2 - A-'Qi-Gd^i^g^'^'^)) (6.48) 



50 



and the approximate function 

u{t; to){z) = ip{z - ziito); c) + 5vi{t; to){z - 2:1 (to)) + Sg{z - Z2{to)), (6.49) 

where 

<f/l,t/l> ' <U2,U2> 



Applying the RG equation to u(t; to){z), one has 

+ 5 \{-C - a}ui'^ + {C- P}ui'^] + 0(6^) 



du 



to=t 



(6.51) 



where u[^\z) = Ui{z — -2i(to)) = dz<f{z — Zi{to)) and U2''' is similarly given. Since 



(1) 



f/« = d,ip'-^\ (|63l|) leads to 



Recalling that 



and 



6a = a 



513 = 13 



= -5C{t) -Sa, C = (3. 

< tJi, -6dA^Sg^^^) > _ <Uu -6gy ((^v?W) > 

<Ui,Ui> <Ui,Ui> 

< U2, -QdA^^g^^^) > < U2, -69,,(^(/pW) > 



(6.52) 



<U2,U2> <U2,U2> 

where Lp^^\z') = ip{z' — h), and putting 5C{t) = Ci{t), we finally obtain the equations of 
motion governing the position of the kink and the speed of it, 

<L^i,-6g,.(v.y.W)> . < U2, -GdAvV^'^) > 
Zi = — Ci H ~ , Ci = . (6.5-3) 



<u,,u,> 



< U2,U2 > 



Similarly, one can readily obtain the equation for £2. This is the main result in this 
subsection. 

Eliminating Ci from the above equations, one can compare our result with those given 
by Ei and OhtaHH]. We first notice that 



d ^ d <[/i,-6g,,((pyW)> 2^ 

—a = = (J{d ' 

dt dt < Ui, Ui > 



(6.54) 



because a depends only on h = Z2 — Zi and h = 0{6). Then we have 



Zl 



<^„6g.,(^V.W)> ^ _ig^5/2^-^(..-.o + 0(6% 

< U2,U2 > 



(6.55) 
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and similarly 



Z2 



(6.56) 



which coincide with the result by Ei and Ohta|^. In short, we have derived the equation 



describing the soliton-soliton dynamics in the KdV equation given by Ei and Ohta with 
fewer ansatz on the basis of the RG method. 



As noticed above, the system (|6.53|) have more information for the soliton-soliton 
dynamics than (|6.55|) (or (|6.56|) ) because (|6.53|) gives the the change of the velocity C\ 
as well as the speed of the position z\. It is intriguing to see an explicit form of (|6.53| ), 
which will provide more detailed information of the dynamics of the solitons. 



The right hand side of the equation for C\ in (|6.53|) is already given in ( |6.55|) . Let us 
calculate the r.h.s. of Z\. The adjoint eigenfunction Ux is given by 

Ui = dcf{s; c)ds + 6*^3, 

J — oo 

where 6 = -M1/M2 with Mi = 1/2 ■ (^J dc^pds^ and M2 =< (p, dc^p Note that 

(C/i,C/2) = holds. M2 is evaluated to be ^/c/2^. We here evaluate Mi. We first 
notice that dc^p = ip/c + zdz(p/2c. Then by a partial derivative, we have 



/ dcfds = 77- f / (fids + 



Zip 



I roo I 

hence Mi = l/2c and 

9 = -l/c^c. (6.57) 



Since Ui is related with f/2(= ^) through (|6.42|) , we see that 

Qifd.Ui = cd.Ui - d^Ui + ip. 

Then performing a partial derivative, we have 

(f7i,-69,(<^<^('^))) = (c9,f/i-9,3f/i + <^,v9W) 

~ 2ce~^'^ / e^^{ca,f/i - 9^f/i + ip)dz, (6.58) 

J ~oo 

where we have used the relation ip^^^{z) ~ 2cev^(^-'^) for h >> 1. The remaining integral 
can be performed using the the asymptotic forms of the derivatives of ?7i, 

dJJx ~ (4 - v^;z)e-^^ dllJx - -v^(5 - v^2)e"^". (6.59) 
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Thus, for L >> 1 we have 

e^^'d.Uidz = e^'^Ui{L) ~ y/^ e^^'Uidz, (6.60) 

J ~oo J~oo 

e-^^'&lUidz = -v^(9-2v^L) + ce^^f/i(L) -cv^ r e^^t/idz. (6.61) 

J —oo J — oo 



Finally, the third term of ( |6.58| ) is evaluated to be 



e^^'ifidz = e^'sech^ zdz {L' = VcL/2), 

J — oo J — oo 



2v/^(2L' + log(l + e-^') + ;^- 
2Vc(2L' - 1) = 2v^(v^L - 1). (6.62) 



Inserting (lOoD ~ (|;6|) into ( [g^SSD , we have 



L^oo J-oo 

= 14v^. (6.63) 
The normalization integral is evaluated to be 

( f/i, f/i ) = - ( / dc(pdz, d,ip) = - { dc^p, V ) = -v^/2. 

J—oo 

Thus we finally obtain 

<[/..-6a,(^^(^))> ^ 

< Ui,Ui > 



and hence the equations for zi and Ci is 



ii = -Ci - 28ce-^(^2-^i)^ 
Ci = I6c5/2e~^(^2-^i)^ 



(6.65) 



together with a similar equation for Z2- Notice that Eq. ( |6.65| ) is consistent with the order 
estimate of Eq. (|6.55| ) . 
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7 Brief summary and concluding remarks 



We have formulated the RG method as a powerful tool for reduction of evolution equa- 
tions in terms of the notion of invariant manifolds, starting from the exact Wilson RG 
equation. We have given an argument as to why to should be set t = in the perturba- 
tive RG method. We have shown that the perturbative RG method constructs invariant 
manifolds successively as the initial value of evolution equations; the integral constants 
in the unperturbative solution constitutes natural coordinates of the invariant manifold 
when the linear operator A in the evolution equation has no Jordan cell. When A has 
a Jordan cell, there is a slight complication because the dimension of the invariant man- 
ifold is to change by the perturbation. The RG equation determines the slow motion 
of the integral constants in the unperturbative solution on the invariant manifold. We 
have worked out several examples to demonstrate our formulation. We have emphasized 
that the underlying structure of the reduction by the RG method completely fits to the 
universal structure elucidated by Kuramoto[|l^ a decade ago. The prescription suggested 



by the present formulation has turned out to be the same as that adopted in [p5| , |26[| . We 
have applied the method to interface dynamics such as kink-anti-kink and soliton-soliton 
interactions in the latter of which a linear operator having a Jordan-cell structure appears. 

In the present work, actual calculations are all based on the perturbation theory, 
although we have started the formulation with the exact flow equation. Recently, variants 
of the exact RG equations or flow equations are applied to various problems such as the 
chiral symmetry breaking in QCD[0, the Bose- Einstein condensation in Alkali atoms 
|6l|| and so on. There, in stead of the perturbation theory, a truncation of the functional 



spaces are usually employed as a practical method of calculations. One may thus imagine 
that a possible practical method for applying the non-perturbative RG method could be 
a scheme similar to, say, the Galerkin method] 6^. Such a non-perturbative RG method 



should be most interesting for analyses of partial differential equations. An attempt to 
apply a kind of non-perturbative RG equation to partial differential equations is given in 



We have tried to formulate the RG method so that the mathematical structure of 
the method becomes as transparent as possible. Then the present work could be a basis 
for clarifying a possible relation between the RG method and another powerful theory 
for reduction of evolution equations called Whitham's averaging method |]6^: The latter 
method has been successfully used to extract the equations describing modulations of 
dispersive non-linear waves: Modulations of the phase function are given by slow variables 
which are governed by so-called Whitham equations. One may thus imagine that the 
modulations as described by Whitham equations may also be given by the RG method as 
it gives the amplitude and phase equations. Or more strongly, Whitham equations might 
be derived as RG ones||65||. Furthermore, Whitham gave a foundation of his equations on 



the basis of a variational principle for the actions, which is reminiscent of the fact that 
Wilsonian RG is formulated for effective actions which admit a variational principle. 
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Appendix A An efficient operator method of solution suitable 
for the RG method 



In this Appendix, we shall summarize rules for obtaining appropriate special solutions 
of non-homogeneous equations appearing in the higher orders in the perturbative RG 
method. The essential point of the present method is to consider the equations with 
initial conditions at t — \/to; the initial values are determined so that the terms in the 
special solutions disappear which either could be " renormalized away" by a redefinition 
of the integral constants in the unperturbed solution or are the ones describing a rapid 
motion. We shall show that the simple rules to write down the special solutions with such 
initial conditions can be summarized as an operator method. 

Let us consider the special solution to the equation given by 

{dt-A)u{t;to) = F{t), (A.l) 
with the initial condition ai, t — to, 

u{to,to) = W{to). (A.2) 

The solution reads 

u{t; to) = e^(*-*°) W^(to) + e^* T dse-^^F(s). (A.3) 
Let ^4 be a semi-simple matrix with eigenvalues Aq, {a — 1, 2, ...); 

AUa = XaUa. (A.4) 

The case where A has a Jordan cell will be considered later. 
Let 

F{t) = e^«*t7^ + e^^*C/« + e^«*C7^, (A.5) 
with Xa ^ Xp, then the solution is evaluated to be 

u{t;to) = e^(*-*°) 



Ap — A Aq — A 



+ (t - to)e""*t/„ + -^e^/'*!/, + —^e^^^Up. (A.6) 

Ap — A Aoc — A 

The first line suggests that the initial value should be chosen as 

^(^o) = Y^^^''"Uo^ + —^e^-'^Up, (A.7) 

Ap — A Aa — A 
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where the first term corresponds to the one which could be renormalized away by a 
redefinition of the unperturbed solution and the second term to a rapid motion. Thus we 
end up with the special solution given by 

u{t; to) = {t- to)e^-'U^ + -^—e^'^'U^ + -^—re^-'Uf^. (A.8) 

Ap — A Aa — A 



Then one may summarize the results as rules of an operator method for obtaining 
special solutions as follows; 

^ e^'Ua, (X^Xa), (A.9) 



A — Aq, 



dt-A " dt-X,^ 

= {t-to)e^-'U^. (A.IO) 

Similary, one can verify that 

Ot- A Of - Xa 

^ (t - to)"+ V«*C7,. (A.ll) 



n+1 



Furthermore, when A 7^ Aq,, 



dt-A' " dr-A 

1 



Ato fin \t t t 



T=t—to 



T=t-to 



(A.12) 



Hence, for example. 



^ -{t-to)e^'Ua = -^{{t-to)-^^}e^'U^, (A.13) 



P_^v" A-^^^ X-A- 



dt-A' -^^'--^ \-A^' X-A' {X-Af 

where A may be replaced with A^. 

Next, we consider the case where A has a semi-simple zero eigenvalue; 

AUq = 0. (A.15) 

Let P and Q be the projection operator onto the space spanned by Uq and its orthogonal 
compliment, respectively. If G is a constant vector, then one can easily verify that 

Si^^'^ = jPG = jjsPG=it-t,)PG, (A.16) 

^QG = ^QG. (A.17) 
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Similarly, 



df-A 



1 OO OO 1 

f{t)QG = —Y^(A-%rQG=-Y,f-\t)—QG, 



■'t „=o n=0 

with b eing the n-th derivative of /(t). Thus, for example, 

^ {i-i^fPG = ^(t-to)-+ipG, 



dt-A 



n+ 1 



1 



(t-to)QG 



(t - t^fQG 



1 



1 



A^ A3 



(A.18) 
(A.19) 

(A.20) 
(A.21) 
(A.22) 



Finally we consider the case where A has a two dimensional Jordan cell; 

AUx = 0, AU2 = C/i. (A.23) 

The adjoint has also a Jordan cell; 

A^XJx = 0, A^U2 = ill. (A.24) 

The adjoint operator A^ is defined by {V,AU) = {A^V,U), where {V,U) is the Her- 
mitian inner product. We define the projection operators P and Q onto the subspace 
{U 1,112} and its orthogonal compliment, respectively. We suppose that the following 
normalization condition is satisfied; 

(C/2,j7i) = l, (C/i,t/2) = l. (A.25) 

A vector U in the P-space is decomposed as U — {U2,U)Ui + {Ui,U)U2- 

Let G be a constant vector, then one has. 



1 



-f{t)PG 



j:{dr'ArmPG, 



t n=0 
-j^ 00 



^{drf{t))A-[{U2, G)U, + (C7i, G)U2], 



jm{{U2, G)U, + (t/i, G)U2} + jjmiJ^, G)U,], 
'dsf{s)PG+ f ds r ds'f{s'){Ui,G)Ui. 



Thus, for example. 



to 



PG 



to ''to 



(A.26) 



{t-to)PG+^{t-toy{Ui,G)U,, 



d,-A 



{t - toTPG 



1 



{t - toT^'PG, 



(n + 1) 



{t-tor+^Ui,G)Ui. 



(n + 2)(n + l) 



(A.27) 



(A.28) 
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The formulae involving QG are the same as those in the semi-simple case. 

The extension to the case where A has a higher dimensional Jordan cell is easy. For 
instance, when A has a three-dimensional Jordan cell such as 

AC/i = 0, AU2 = Ui, AU3 = U2, (A.29) 

one can easily verify that 

-\f{t)PG = f dsf{s)PG + /■* ds r ds'f{s'){{U2, G)U, + (C/i, G)t/2}, 

Ot — A Jto Jto Jto 

+ f ds f dsi r ds2f{s2){Ui,G)Ui, (A.30) 

where the adjoints satisfy 

A^ili = 0, A^U2 = ill, A^Us = U2. (A.31) 
The normahzation condition reads (t/3, Ui) — {TJ2, U2) — (t/i, C/3) = 1. 
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Appendix B An elementary method to derive the approxi- 
mate solution for the double-well potential 



The first integral of tlie Newton equation with the initial condition x(0) = for the 
double-well potential reads by 

rx/V2E dy 

Expanding the integral in a Taylor series, one readily obtains 

t = (1 _ ^eE)Smh-'X + e- ^' + ^^ , (B.2) 

with X = x/V2E up to 0{e^). Here Sinh'^X = In |X + y/X^ + 1|. The point of the 
present method is to notice that the above equation is equivalent to the following equation 
up to O(e^); 

« = Sinh-iX + ^^l±^, (B.3) 

with M = (1 + ^eE)t. One may solve ( |B.3| ) perturbatively and obtains 

3 . , e 



x{t) = V2E{1 - -eE)smhu - -{V2E)^ sinh^ u. (B.4) 
4 8 

Putting -\/2E{(l — ^eE) = C, one ends up with 

x{t) = C sinh at - ^C'^ sinh^ at, (B.5) 

with a = 1 + feC^ up to O(e^), which coincides with the result given in subsection 3.3 in 
the text. 
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Appendix C The period of the Lotka-Volterra equation 



The Lotka-Volterra equation admit periodic solutions. An approximate but globally 
valid periodic solution was explicitly constructed by the RG method by one of the present 
author, [p^ The main purpose of this Appendix is to show that the solution constructed 
there gives the period which coincides with that obtained by Frame [p5|] in a quite different 



approach. To make the argument self-contained, however, we shall repeat the RG analysis 



but in a mathematically more simple way than that given in . 

Introducing the new variable u = *(,^, rj) by 

x={b + eO/e', y = a/e + T], (C.l) 
Eq. (|5.19| ) is reduced to the following one: 

|-Lo).^-.e.(_\), (C.2) 



^={i)^ ^o=f° n). (C.3) 



7] J' \a 



where 



Lq has the eigenvalues A = ±iV ab = ±iuj, with the corresponding eigenvectors given by 
Ui = —ioj/h) and U2 = respectively. Here A* denotes the complex conjugate of 
A. To make the following calculation as transparent as possible, we first transform the 
equation to the form where Lq is diagonalized. Then one finds that the vector equation 
( |C2| ) is reduced to a scalar equation 

(!L.^..^(,.,, ,C.4) 

where a = 1/2 • (1 — ih/oj) and z(t) = 1/2 ■ {(^ + ib/u ■ rj). This equation apparently has a 
much simpler form than that treated in |^ . 



Now we try to solve (|C.4|) around t ~ Vto by the perturbation theory by expanding 
z = zq + ezi + e^Z2 + o(e^), with the initial condition zit^to) = W(to). Wito) is also 
expanded as W{to) = Wo{to) + eWi{to) + eW2(to) + o{e^). 

A simple manipulation gives the solution as follows; 

¥ \ 3 
o lap f 1 2iuj 
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|_ 1^3g-3..i + f^CC\t - to)e^"*} . (C.5) 
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Now the RG equation -§^it;to) 1^^=^= gives 

Here we have used that |ap = i(l + ^). The general solution z{t, to) is now given as the 
initial value W(t) by construction z{t) = W(t). 

Since |C(t)p = const, as easily verified from ( |C(j| ), one may put C = ^e*^, with A 
being a real constant. Then ( |C.6|) implies that 6 = — ■ (1 + b'^/uj'^)uj, hence 

e{t) = -£2^2/2462 . (1 + byuj^)ujt + e. 

For a calculational convenience, let us define Q hj <d{t) = ut + 6 with 



which implies that C (i(:)e"^* = A/2 - (sin Q — i cos G). Then we have finally the components 

'KW.'K*)), 



Although these expressions are seemingly different from those given in p6[ , they coincides 



with each other up to o(J^\ In fact, by the redefinition of the constant variables ^4(1 
e^^^^) ^ A and ^ - e^^^ ^ ^, ([CJ) and (^) reproduce the result in |2|, up to 



It is remarkable that we have obtained the modified angular velocity cD depending on 
the amplitude. A long ago. Frame gave the period of the motion. Here, we shall show 
that the period Trg = 27r/cl; coincides with that given by Frame up to o{e^). 

Our period reads 

^ 27r 27r , A 1,, , ^ 

Trg = — , , 2 , = ^rS^ + ^ - + 7; > + ' c-^o 



where use has been made of a; = v ab. On the other hand, the period Tpr given by Frame 
reads in his notations 

- + TUT + ^ + ■ ■ ^^-^^^ 
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Here, ci, c are identified with the variables expressed by our basic variables as follows 
9 a + h 



24 



(C.12) 



We remark that the correspondence between the variables of ours and Frame's is summa- 
rized as follows: 

a ; h ^ ; e ; e' 

6 a 
a; ^ A?"! ; y ^ N2 ; — ^ rii ; - ^ n2. 

e t 

It should be remarked here that (? is actually a constant because our system has a 
conserved quantity 61nx — e'x + alny — = const. In fact, one can easily verify that 

Thus, 

Tpr = + c^c^} + 0(6^) = + ^4(? + -)} + (C.15) 



/a6 Va6 24 ab b a 

From (^) and (^), one sees that ^ = + o(l), 77=^^^ + o(l). Hence 



27r 1 1 4^ , 2x 

= Trg + o(£2). (C.16) 
This is what we wanted to show. 
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